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Nichols,  Todd  W.  (Ph.D.,  Electrical  Engineering) 

Measurements  and  Calculations  of  Loss  Due  to  Finite  Conductivity  and  Surface 
Roughness  in  Microstrip  Lines 
Thesis  directed  by  Assistant  Professor  John  M.  Dunn 

Finite  conductivity  and  surface  roughness,  two  often  significant  causes 
of  loss  in  microstrip  lines,  are  examined  in  this  dissertation.  First,  actual  data 
are  gathered.  Gold  lines  are  fabricated  on  semi-insulating  gallium  arsenide  sub¬ 
strates.  Losses  are  measured  for  the  three  canonical  types  of  surface  roughness 
in  rectangular  cross  section  lines:  smooth  lines,  lines  with  grooves  parallel  to  the 
line  length,  and  lines  with  grooves  transverse  to  the  line  length.  Multi-line  thru- 
reflect-line  calibrations  are  performed,  and  on- wafer  measurements  are  made  on 
both  long  lines  and  two-port  gap-coupled  resonators.  Data  are  compared  with 
other  methods  of  loss  calculation,  including  commercial  software. 

The  second  major  effort  is  solution  of  an  exact  integral  equation  for 
the  current  distribution  in  the  strip  conductor  of  a  microstrip  line  on  a  per¬ 
fect  dielectric  above  a  perfect  ground  plane.  The  strip  conductor  may  have 
finite  conductivity  and  arbitrary  edge  shape.  A  quasi-transverse  electromag¬ 
netic  mode  is  assumed  and  the  moment  method  is  applied  to  a  discretized 
strip.  The  resulting  matrix  equation  is  a  nonlinear  complex  eigenvalue  problem 
for  which  the  eigenvalues  are  the  propagation  constants  of  the  possible  modes. 
The  real  part  of  the  dominant  eigenvalue  is  the  desired  loss.  The  case  with  air 
dielectric  (e^  =  1)  is  investigated. 
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CHAPTER  I 


INTRODUCTION 


1.1  Motivation 

The  lowly  microstrip  transmission  line  has  become  firmly  entrenched 
in  such  applications  as  hybrid  microwave  integrated  circuits  (MIC  s)  on  ceramic 
and  polytetraflouroethylene  (PTFE.  or  Teflon)  substrates,  and  monolithic  mi¬ 
crowave  integrated  circuits  (MMIC's)  on  silicon  and  gallium  arsenide  (GaAs). 
We  live  in  an  imperfect  world,  and  microstrip  is  also  subject  to  imperfection- 
the  imperfection  of  loss:  loss  due  to  radiation,  imperfect  dielectrics,  and  im¬ 
perfect  conductors.  Radiation  losses  arise  from  discontinuities,  dielectric  losses 
are  due  to  finite  resistivity  or  leakage  of  the  substrate,  and  conductor  losses- 
the  subject  of  this  dissertation-are  the  result  of  finite  conductivity  and  surface 
roughness  of  the  line's  conductors. 

Figure  1.1  shows  the  geometry  of  a  “perfect"  microstrip  with  perfect, 
lossless  dielectric  and  conductors  and  zero-thickness  strip  conductor.  The  first 
step  toward  reality  is  to  acknowledge  that  the  strip  has  some  nonzero  thickness 
t.  The  substrate  has  a  loss  tangent  -which,  for  good  microwave  substrates,  can 
usually  be  ignored  (otherwise,  it  is  not  a  very  good  substrate).  The  next  step 
is  to  include  effects  of  a  finite  conductivity  a.  In  this  dissertation,  the  conduc¬ 
tor  losses  dominate,  both  in  the  substrates  considered  and  in  the  phenomena 
studied;  therefore  the  loss  tangent  will  be  neglected.  Finally,  no  metallic  or 
ceramic  surface  is  perfectly  smooth.  So  we  introduce  a  degree  of  surface  rough¬ 
ness  (usually  given  as  a  root-mean-square  amplitude  of  the  cross-section  of  the 
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perfect  zero-thickness  conductor 

W  — ^ 

perfect  dielectric 
perfect  conductor 

Figure  1.1  "Perfect"'  microstrip  model 

substrate  height  variation),  and  arrive  at  the  "real-world*"  microstrip  depicted 
in  Figure  1.2. 

Conductor  and  surface  roughness  losses  in  microstrip  have  been  exam¬ 
ined  in  great  detail  previously  (see.  for  example,  [1-7]);  Wheeler  s  incremented 
inductance  rule  has  especially  seen  considerable  use  [1].  This  document,  though, 
presents  measurements  made  on  a  controlled  rough  surface,  as  well  as  exami¬ 
nation  of  an  integral  equation  which  is  not  usually  used  for  loss  calculations. 

1.2  Effects  of  Finite  Conductivity  in  Microstrip 

The  conductor  loss  in  a  transmission  line  supporting  a  cjuasi-TEM 
mode  is  usually  calculated  with  a  perturbational  technique  [S-9].  Using  trans¬ 
mission  line  theory,  we  can  write  a  power-current  definition  of  the  line"s  char¬ 
acteristic  impedance: 


P  =  Zo\I\- 


(1.1) 
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imperfect  conductor 


rough  interface 


Figure  1.2  Microstrip  with  real-life  imperfections 


where  Zo  is  the  characteristic  impedance  and  I  is  the  current  flowing  down 
the  line.  Substituting  in  the  I+e~°-  form  of  a  wave  traveling  in  the  positive  :: 
direction,  we  have 

P=ZoiJ+pe-2°^=Poe-2-  (1.2) 

where  Po  is  the  power  at  r  =  0  and  is  the  direction  of  the  line.  Now  the  rate 
of  power  loss  per  unit  length  at  some  position  r  is  given  by 

PL  =  -^=2aPoe-2-;  (1.3) 

oz 


therefore 


nepers/unit  length 


(1.4) 


This  is  the  classic  relation  from  which  most  derivations  start. 

The  next  step  is  to  find  and  evaluate  expressions  from  the  two  powers 
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in  Equation  (1.4):  usually  a  Poynting  vector  approach  is  taken: 

P  =  [  E  X  H* -a^dS  (1.5) 

2  Js 

and 


where  Ci  and  C'2  are  the  contours  of  the  two  conductors  and  S  is  the  line's  cross 
section.  P  is  the  power  flowing  down  the  line.  Pl  is  the  (lost  or  dissipated) 
power  flowing  into  the  conductor.  These  two  equations  are  still  general. 

Conventional  approaches  simplify  (1.5)  and  (1.6)  with  the  usual  as¬ 
sumption  that  the  Leontovitch  impedance  boundary  condition  is  valid  (that 
is,  the  line  thickness  t  is  much  greater  than  a  skin  depth,  and  the  conductor 
surfaces  are  smooth).  That  means  that  the  fields  internal  to  the  conductor 
do  not  need  to  be  calculated.  Unfortunately,  for  modern  microstrip  lines  this 
approximation  is  often  not  valid.  Other  more  drastic  measures  are  called  for. 

1.3  Effects  of  Surface  Roughness  in  Microstrip 

It  is  accepted  that,  since  metals  are  known  to  have  finite  conductivity 
(for  gold,  a  =  4.1E7  S/m  and  for  copper,  a  —  5. SET  S/m),  there  will  be  some 
impediment  to  wave  transmission  and  that  it  must  be  calculated.  This  is  usually 
taken  into  account  by  relating  power  dissipation  to  the  current  in  a  conductor 
through  Ohm’s  Law  {E  =  crj).  This  means  that  the  integrals  in  (1.5)  and 
(1.6)  have  to  be  evaluated,  as  previously  mentioned.  But  what  about  surface 
roughness?  Surely  we  can  see  that  the  substrates  (and  metals)  we  are  using  are 
smooth?  Let  us  examine  them  somewhat  closer.  As  we  do,  let  us  bear  in  mind 
that  the  skin  depth  8  (in  gold)  can  range  from  about  0.97//  at  6  GHz  to  about 
0.38//  at  40  GHz  (the  frequency  range  examined  in  this  work).  Remember,  this 
is  the  dimension  which  determines  whether  or  not  a  surface  really  is  rough. 
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Only  two  substances  (substrate  and  conductor)  make  up  microstrip, 
so  only  those  two  substances  can  be  rough.  First  we  look  at  substrates.  Mi¬ 
crostrip  lines  are  fabricated  on  three  main  classes  of  dielectrics:  PTFE,  which 
is  easily  worked  and  inexpensive;  ceramics  (predominately  aluminum  oxide), 
which  are  hard  to  work  and  fairly  expensive;  and  semiconductors,  which  are 
used  in  MMIC’s  and  may  have  a  wide  cost  and  workability  range.  PTFE's  sur¬ 
face  roughness  depends  on  processing  and  the  amount  and  type  of  ceramic  fibers 
mixed  in  with  it,  and  may  range  from  0.3;u  to  1.0/.t  [10].  For  ceramics  we  will 
concentrate  on  alumina  (AI2O3)  since  it  is  the  most  common  [11].  Alumina  is 
available  in  several  purities,  and  the  purity  affects  post-firing  surface  roughness. 
Of  course  the  surface  can  be  polished,  but  since  alumina  is  9  on  the  Mohs  scale 
of  hardness,  it  is  expensive  to  do  so.  Fired  unpolished  99.5%  alumina  may  have 
a  surface  roughness  from  l/J.  to  as  high  as  S/U,  depending  on  forming  technique 
[12-13],  whereas  96%  alumina  (cheaper)  may  have  a  gross  surface  roughess  of 
20p.  Polishing  can  reduce  this  roughness  to  0.04/u  on  99.5%  alumina  [14].  So 
both  PTFE  and  unpolished  ceramic  substrates  can  lead  to  real  problems.  Semi¬ 
conductors  are  usually  supplied  from  the  wafer  manufacturer  optically  fiat,  so 
there  is  no  substrate  roughness  problem  there. 

Now  we  will  consider  the  metal.  Conductor  can  join  substrate,  in 
which  case  the  rougher  and  tougher  material’s  roughness  dominates,  or  it  can 
join  air  (the  top  of  the  strip)  -where  the  metal’s  roughness  (process  dominant) 
matters  more.  Copper  is  usually  rolled  or  electrodeposited;  rolled  copper  surface 
roughness  may  range  from  0.1/u  on  the  substrate  side  to  1.0//  on  the  air  side 
[15].  Electrodeposited  copper,  on  the  other  hand,  is  a  factor  of  ten  rougher, 
which  would  seem  to  preclude  its  use;  however,  it  adheres  much  better  to  the 
substrate  [16].  Gold  is  also  sometimes  used  on  ceramic  substrates.  Gold  is 
also  the  primary  metal  used  on  semiconductors.  It  is  either  evaporated  or 
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electroplated  onto  the  substrate;  the  electroplating  usually  takes  place  over  an 
evaporated  seed  layer  that  may  be  a  few  thousand  Angstroms  thick.  Gold 
electroplated  onto  wafers  made  at  Wright-Patterson  Air  Force  Base  was  seen 
to  have  a  surface  roughness  of  about  0.1  jU  to  0.4//  as  determined  from  Dektak 
measurements. 

One  way  to  approximate  surface  roughness  to  make  it  amenable  to 
analysis  is  to  look  at  grooves,  both  parallel  to  and  transverse  to  current  flow 
in  the  strip.  Morgan  [17]  made  some  early  studies  on  the  effects  of  grooves  on 
attenuation.  For  parallel  grooves  rvith  width  of  channel  ecpal  to  width  of  mesa, 
as  the  groove  depth  to  skin  depth  ratio  increases,  the  attenuation  increases  also. 
Conceptually,  this  is  because  the  current  is  trying  to  flow  through  a  layer  which 
is  reduced  in  area  (by  the  grooves),  so  there  is  less  metal  to  flow  through.  For 
transverse  grooves,  as  this  ratio  increases,  the  path  length  for  the  current  flow 
effectively  increases,  so  the  attenuation  goes  up.  If  groove  depth  equals  groove 
width,  the  loss  could  double.  Of  course,  in  real  life  the  surface  roughness  is  not 
this  ideal,  but  it  is  obvious  that  surface  roughness  -  especiallj'  that  on  the  order 
of  a  skin  depth  -  can  noticeable  increase  the  loss. 

1.4  Arrangement  of  the  Dissertation 

The  first  major  task  in  this  research  was  the  design  and  fabrication 
of  a  gallium  arsenide  (GaAs)  MMIC  with  gold  conductors.  Chapter  2  details 
the  items  contained  on  the  wafer  and  the  techniques  used  to  process  it.  This 
effort  was  spent  mostly  in  the  clean  room;  there  is  little  to  show  for  it  but  the 
wafers  themselves  -  and  a  fairly  well-debugged  process  flow.  Producing  these 
wafers,  however,  consumed  much  time.  Trial  and  error  associated  with  the 
initial  wafers  (as  detailed  in  Appendix  A)  was  a  factor,  as  well  as  the  fact  that 
the  final  wafer  design  possessed  via  holes  and  required  backside  processing. 

Chapter  3  investigates  the  loss  measurement  methodologies  employed 


to  actually  acquire  data.  Appendix  B  describes  the  hookup  and  use  of  the 
measurement  equipment  itself,  since  this  facility  —  which  is  capable  of  probing 
on-wafer  active  devices  to  40  GHz  —  took  form  mostly  through  this  research; 
several  good  measurement  habits  were  developed  and  documented.  Since  net¬ 
work  analyzer  calibration  is  an  integral  part  of  these  procedures,  calibration  - 
particularly  thru-reflect-line  (TRL)  calibration  -  is  described.  Measuring  the 
loss  of  long  lines  directh'  and  deducing  the  loss  from  resonator  parameters  are 
the  two  popular  methods  for  loss  measurement,  and  these  two  methods  are 
scrutinized.  Finally,  the  ground  plane  loss  must  be  calculated  and  accounted 
for.  since  the  ground  plane  is  not  roughened  in  this  experiment. 

The  hard-won  data  are  presented  and  discussed  in  Chapter  4.  After  a 
thorough  discourse  upon  this  data,  results  from  traditional  approximate  formu¬ 
las,  other  researchers,  and  commercial  software  calculations  are  compared  with 
these  numbers. 

Chapter  5  marks  a  dramatic  shift  in  subject  matter.  The  solution  of  an 
exact  integral  equation  for  the  current  distribution  in  a  finitely  conducting  strip 
of  arbitrary'  cross  section  above  a  perfect  dielectric  upon  a  perfect  ground  plane 
is  outlined.  Appl3'ing  the  method  of  moments  results  in  a  complex  nonlinear 
eigenvalue  problem  in  which  7  =  a  is  the  eigenvalue.  Starting  values  of 
7  are  provided  bj'  traditional  approximate  formulas.  The  special  case  of  air 
dielectric  (cr  =  1)  is  first  investigated;  this  approach  eases  and  modularizes  the 
solution.  First,  this  case  involves  simple  contour  integrals  instead  of  the  more 
difficult  to  handle  Sommerfeld  integrals.  Second,  once  the  easier  air  dielectric 
case  is  working  properly,  the  dielectric  can  be  incorporated  by  changing  the 
Green’s  function  used.  Finallj',  since  the  air  dielectric  case  is  easier  to  solve,  it 
is  easier  to  get  working  initially.  With  some  forethought,  the  framework  easily 
accommodates  the  nonunity  dielectric  constant  later.  There  was  not  enough 
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time  to  complete  addition  of  the  dielectric,  so  only  the  air  dielectric  case  is 
presented  in  detail. 


CHAPTER  II 


DESIGN  AND  FABRICATION  OF  A  WAFER 
SUITABLE  FOR  LOSS  MEASUREMENT 

How  may  one  measure  the  loss  due  to  finite  conductiAnty  and  sur¬ 
face  roughness  in  microstrip?  Obviously  the  question  necessitates  the  use  of 
some  form  of  microstrip  transmission  line.  ith  that  decided,  the  technology 
(substrate  type,  metallization  type,  and  fabrication  method)  must  be  selected. 
But  yet  another  question  intrudes:  How  does  one  make  a  rough  surface?  This 
answer  will  surel}"  influence  the  decision{s)  about  technology. 

2.1  Design  goals 

Simply  stated,  the  goal  was  to  measure  the  loss  in  microstrip  due  to 
finite  conductivitj’  and  surface  roughness,  in  the  frequency  range  where  the 
strip  thickness  and  surface  roughness  are  both  on  the  order  of  a  skin  depth. 
The  effects  of  no  surface  roughness  (on  smooth  lines),  grooves  running  parallel 
to  the  line  length  (parallel  grooves),  and  grooves  running  transverse  to  the  line 
length  (transverse  grooves)  were  to  be  studied.  These  are  the  three  different 
canonical  problems.  Any  given  surface  roughness  can  be  expressed  as  a  linear 
combination  of  parallel  and  transverse  roughness. 

Under  the  terms  of  the  Air  Force  fellowship,  the  clean  room  facilities  at 
the  Avionics  Laboratory  at  Wright  Research  and  Development  Center  (WRDC) 
located  at  'V^^ight-Patterson  Air  Force  Base  in  Dayton,  Ohio  were  available. 
With  this  availability  and  with  no  readily  discernable  method  for  making  a 
carefully  controlled  rough  surface  in  more  tradional  microwave  substrates  (i.e.. 
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alumina  and  Teflon),  the  decision  was  made  to  fabricate  gold  lines  on  semi- 
insulating  gallium  arsenide  (GaAs)  wafers. 

As  noted  earlier,  surface  finish  and  line  thickness  are  of  concern  when 
their  characteristic  lengths  are  on  the  order  of  a  skin  depth.  The  skin  depth 
S  for  gold  at  40  GHz  is  about  0.377  n  and  increases  to  0.972  as  frequency 
decreases  to  6  GHz.  If  the  grooves  are  about  one  to  one-half  micron  deep,  then 
a  strip  conductor  2.5  /<  thick  will  both  be  thicker  than  twice  the  groove  depth 
(so  that  there  will  always  be  metal  across  the  width  of  the  line)  and  be  thin 
enough  to  be  only  a  few  skin  depths  thick. 

Fifty  ohms  is  the  most  common  system  characteristic  impedance,  and 
100  n  is  the  most  common  GaAs  MMIC  substrate  thickness.  Also  the  substrate 
must  be  thin  in  order  to  suppress  undesired  modes  at  frequencies  to  40  GHz. 
Thus  if  a  fifty-ohm  line  is  desired  on  100  n  GaAs  (whose  relative  permittivity 
is  12.9),  its  width  must  be  about  73  /z;  this  value  was  ca;lculated  with  standard 
design  formulas  [3,  5,  12].  For  a  line  of  these  dimensions,  the  lowest-order  TM 
mode  cutoff  frequency  is  about  292  GHz,  and  the  lowest-order  transverse  strip 
resonance  mode  occurs  at  about  370  GHz;  therefore  this  substrate  is  sufficiently 
thin.  Also  vertical  edges  are  desired;  therefore  liftoff  was  the  process  used  for 
fabricating  the  metal  lines.  Likewise  vertical  edges  were  desired  for  the  grooves; 
reactive  ion  etching  (RIE)  for  this  shallow  depth  provides  vertical  side  walls. 
As  a  backup,  all  grooves  on  the  wafer  run  in  the  same  direction,  so  that  if 
chemical  etching  is  ever  desired,  the  anisotropy  of  certain  etchants  [IS]  will  not 
be  a  factor.  That  is,  all  the  grooves  would  always  have  the  same  cross-sectional 
appearance  if  chemical  etching  produced  different  edge  slopes  in  different  crystal 
lattice  directions.  Figure  2.1  shows  the  desired  cross  section  of  a  microstrip  line 
with  grooves  running  parallel  to  the  direction  of  propagation. 

Finally,  there  is  the  question  of  actually  making  the  measurements. 
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Figure  2.1  Cross  section  of  niicrostrip  line  with  parallel 
grooves  (not  to  scale) 

For  less  expensive  and  larger  substrates,  coaxial  connectors  are  soldered  onto 
the  ends  of  the  substrate.  This  approach  is  obviously  infeasible  for  a  wafer 
technology:  on-wafer  prol^ing  is  required.  Some  sort  of  transition  is  therefore 
needed,  because  all  wafer  probes  are  variations  of  coplanar  lines,  ground-signal, 
ground-signal-groimd  (traditional  coplanar  waveguide),  or  other  combinations. 
The  microstrip  ground,  though,  is  on  the  bottom  of  the  substrate  and  is  nor¬ 
mally  inaccessible  from  the  top  of  the  wafer.  Two  approaches  are  available: 
using  pads  on  the  top  surface,  connected  to  the  ground  plane  with  metallized 
via  holes;  or  using  a  large  ground  patch  on  the  top  surface  whose  reactance 
to  gromid  is  minimized  in  some  frequency  band  of  interest  [19].  The  latter 
approach  was  used  first  and  later  abandoned:  it  is  presented  in  more  detail  in 
Appendix  A.  The  resonant  patch  has  the  advantage  of  not  requiring  an  extra 
mask  laver  and  extra  processing  for  via  holes.  However,  it  is  a  resonant  struc- 
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Figure  2.2  Coplanar  waveguide  to  microstrip  transition 

ture  and  contributes  to  radiation  losses.  It  also  requires  much  more  space  on  the 
wafer.  Transitioning  from  wafer  probes  to  microstrip  with  a  via-hole  transition 
is  the  standard  method.  Figure  2.2  shows  the  plan  view  and  cross  section  of  the 
transition  used.  The  pad-to-line  gaps  are  44  /i  as  shown  in  order  to  present  a 
coplanar  waveguide  whose  characteristic  impedance  is  close  to  50fi;  this  value 
was  calculated  with  standard  design  formulas. 

2.2  Wafer  inventory 

■V^Tiat  should  be  included  on  the  wafer?  At  a  minimum,  calibration 
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standards  and  the  structures  to  be  measured  are  required.  There  are  two  pos¬ 
sible  ways  to  make  loss  measurements.  One  method  is  to  use  TR.L  calibration 
with  multiple  lines  (discussed  in  Chapter  3)  to  extract  the  loss;  alternatively,  af¬ 
ter  calibrating  with  TRL,  the  IS21 1  of  a  line  can  be  measured  directly.  The  other 
method  is  to  fabricate  a  resonator  and  use  it  in  one  or  more  resonant  modes: 
by  measuring  the  resonant  frequency,  quality  factor. and  coupling  coefficient, 
one  can  extract  the  loss.  This  method,  though,  only  gives  results  at  discrete 
frequencies.  To  satisfy  adherents  of  both  approaches  and  to  make  comparisons, 
both  methods  are  used  here.  This  means  that  sets  of  two-port  gap-coupled 
resonators  are  present,  as  well  as  sets  of  long  lines.  Also,  since  some  wafer  area 
was  available,  structures  were  included  to  provide  data  for  other  researchers  at 
the  University  of  Colorado.  Data  on  these  structures  are  not  presented  in  this 
dissertation. 

As  a  precaution  against  low  yield,  two  of  every  structure  were  included 
on  the  wafer.  Moreover,  every  line  and  resonator  is  present  in  smooth,  parallel, 
and  transverse  forms.  In  addition,  the  structures  are  spaced  at  least  five,  and 
usually  eight  to  ten,  substrate  thicknesses  apart  to  minimize  coupling  from 
structure  to  structure.  The  design  of  the  particular  calibration  standards  is 
discussed  in  Chapter  3.  Figure  2.3  shows  the  structures  present  on  the  wafer, 
and  Figure  2.4  is  a  view  of  the  entire  wafer.  Table  2.1  is  a  detailed  list  of  the 
structures  present  and  their  dimensions.  Table  2.2  lists  the  frequencies  at  which 
the  resonators  will  give  loss  measurements. 

2.3  Fabrication 

We  will  now  discuss  the  actual  fabrication  procedures  and  process 
steps.  This  design  process  is  essentially  the  same  one  used  in  GaAs  MMIC 
foundries  for  commercial  chips. 
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Figure  2.3  Calibration  standards,  resonators,  and  lines  on 
the  wafer 
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Table  2.1  Complete  wafer  inventory 


TR.L  Calibration  Standards 

4  -  Offset  opens  (500  lu) 

2  -  Zero-length  thru,  T  (1000  //  ) 

2  -  Full- band  line  for  6-40  GHz,  B  (2121  ft  ) 

2  -  Lower-band  line  for  6-13  GHz,  L  (3713 
2  -  Middle-band  line  for  13-27  GHz.  M  (2289  /i) 

2  -  Upper-band  line  for  27-40  GHz,  U  (1769  /u) 

Long  Lines  (2  of  each  line) 

6  -  2.998  cm  Ll_S,  L1_T,  Ll-P  smooth,  transverse  grooves,  parallel  grooves 

6  -  2.498  cm  L2_S,  L2_T,  L2_P  smooth,  transverse  grooves,  parallel  grooves 

6  -  2.098  cm  L3_S,  L3_T,  L3_P  smooth,  transverse  grooves,  parallel  grooves 

2-port  Gap-coupled  Resonators 

6  -  7400  fi  Rl_S,  Rl-T,  Rl-P  smooth,  transverse  grooves,  parallel  grooves 

6  -  6450  //  R2_S,  R2_T,  R2-P  smooth,  transverse  grooves,  parallel  grooves 

6  -  5450  //  R3_S,  R3_T,  R3_P  smooth,  transverse  grooves,  parallel  grooves 

6  -  4650  p  R4_S.  R4_T,  R4_P  smooth,  transverse  grooves,  parallel  grooves 

6  -  3940  fi  R5_S,  R5_T,  R5_P  smooth,  transverse  grooves,  parallel  grooves 

Miscellaneous  Structures 

1  -  13.66  cm  smooth  serpentine  line 

3  -  various  serpentine  lines  with  stubs  (for  other  CU  researchers) 

2  -  double-stub  resonators  (for  other  CU  researchers) 
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Table  2.2  Resonator  frequencies 


Resonator  Number  Length 

1  7400 

2  6450  jj 

3  5450  n 

4  4650  n 

5  3940  // 


Resonant  Frequencies  (GHz) 
7.0,  14.0.  20.9,  28.2,  .34.6 
8.0,  16.0,  23.9,  31.7,  39.5 
9.4,  18.8.  28.1.  37.3 
11.0,  22.0,  32.8 
13.0,  25.8,  38.4 


The  starting  material  was  two-inch-diameter  semi-insulating  GaAs  wa¬ 
fers  approximately  20  mils  thick.  They  were  already  polished  on  one  side  and 
the  major  flat  is  parallel  to  the  <  100  >  crystal  plane.  Fifteen  wafers  were 
available  at  the  beginning.  Table  2.3  gives  the  order  of  the  major  steps  in  the 
fabrication  process,  and  these  steps  will  each  be  elaborated  in  the  rest  of  this 
chapter.  The  processing  steps  were,  unless  noted  otherwise,  carried  out  in  the 
Class  100  ^  clean  room  in  the  Avionics  Laboratory  at  WRDC.  Such  care  is 
taken  with  this  facilitA’’  that  it  is  often  Class  50  or  even  10. 

2.3.2  Grooves  -  first  mask  level 

The  first  step  when  masks  are  in  hand  is  to  fabricate  the  grooA’’es  in 
the  substrate.  Table  2.4  lists  the  steps.  The  general  approach  is  to  put  down 
a  layer  of  photoresist,  expose  the  groove  mask,  develop,  and  use  reactive  ion 
etching  to  form  the  grooves.  Care  must  be  taken  in  exposing  and  developing 
because  the  grooves  are  one  micron  lines  and  spaces.  It  is  very  easy  to  over¬ 
expose  or  underexpose,  and  easy  to  leave  too  much  photoresist  behind  or  to 

^  Class  100  means  that  fewer  than  100  particles  with  diameter  greater 
than  100  are  present  per  cubic  foot  of  air. 
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Table  2.3  Major  processing  steps 


•  groove  lithography 

•  groove  etch 

•  front  metal  lithography 

•  front  metallization 

•  liftoff 

•  lapping  and  polishing 

•  via  hole  lithography 

•  via  hole  etch 

•  back  metallization 

•  demounting 


develop  too  much  away.  Close  inspection  under  magnification  after  developing 
is  mandatory.  Figure  2.5  shows  some  fabricated  grooves.  Note  that  any  parti¬ 
cles  in  the  photoresist  where  grooves  are  to  be  formed  will  distort  the  grooA^e 
pattern  in  that  area.  This  step  is,  therefore,  very  sensitive  to  contaminants  in 
the  photoresist,  in  the  air,  and  in  the  tools  (such  as  pipettes)  used. 

2.3.3  Front  metal  —  second  mask  level 

In  this  step  the  strip  conductors  of  the  microstrip  structures  are  laid 
down.  Table  2.5  gives  the  detailed  process  steps,  while  Figures  2.6  through  2.11 
show  graphically  the  steps  involved.  The  liftoff  process  as  given  here  results 
in  vertical  metal  edges  because  the  base  photoresist  is  thicker  than  the  desired 
metal  thickness,  and  because  the  metal  is  deposited  by  evaporation-an  essen¬ 
tially  line-of-sight  process.  The  overlap  of  the  upper  level  of  photoresist  enables 
solvent  to  be  used  to  dissolve  photoresist  underneath  undesired  metal  so  that 
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Table  2.4  Groove  litliography  and  etching  steps 

Groove  Lithography 

Clean  wafer  with  acetone,  methanol,  isopropyl  30  sec 

Bakeout  at  loO°C  for  10  min 

Spin  on  1400-27  photoresist  at  4000  rpm  30  sec 

Bake  at  100°C  for  5  min  on  hot  plate 

Edge  removal  -  expose  405  nm  UV  10  sec 

Develop  5:1  DI:351  at  500  rpm  30  sec 

DI  rinse  30  sec 

Expose  mask  (grooves)  405  nm  UV  3  sec 
Develop  DI:351  at  500  rpm  13  sec 
DI  rinse  30  sec 

Groove  Etch 

Bake  at  100°C  for  5  min  on  hot  plate  (to  harden  resist) 
O2  plasma  etch  30  sec 
HC1:H20  1:1  rinse  30  sec 
RIE  to  0.8  fi  depth  ±2  fji 


it  will  “lift  off”;  hence  the  name.  Figure  2.12  shows  a  vertical  edge  on  a  2.5  // 
thickness  of  gold  on  a  GaAs  substrate, 

2.3.4  Wafer  thinning 

This  procedure,  though  low  in  number  of  process  steps,  is  dangerous 
-  both  to  the  operator  and  to  the  wafer.  In  this  step  the  wafer  is  reduced  in 
thickness  from  500  fJ.  -  600  //  to  the  100  /u  which  is  the  desired  final  substrate 


Figure  2.5  Scanniug  electron  microscope  photograph  of 
etched  grooves  in  a  GaAs  wafer 

thickness.  This  step  is  dangerous  to  the  wafer  because  GaAs  is  a  very  brittle 
and  fragile  material:  as  it  gets  thinner,  its  fragihty  increases  and  its  toughness 
decreases.  The  substrate  is  also  more  sensitive,  as  it  is  lapped  thinner  and 
thinner,  to  being  scratched  by  an  errant  large  particle.  If  in  the  final  stages  of 
polishing,  when  the  substrate  is  about  105  n  or  so  thick,  a  5  ^  particle  used  in 
the  rough  grinding  gets  in  between  the  substrate  and  the  grinding  wheel,  then 
it  can  scratch  the  substrate  so  deeply  that  smoothing  it  out  will  reduce  the 
substrate  to,  say,  95  /x.  Also  the  grinding  rate  is  not  constant;  it  depends  on 
the  composition  of  the  aluminum  oxide  and  deionized  water  and  Chlorox  slurry. 
It  is  thus  very  easy  to  grind  the  wafer  too  thin.  And  if  the  wafer  is  not  properly 
mounted  to  its  supporting  glass  disk,  parts  of  the  wafer  can  fracture  and  flake 


Table  2.5  Front  metal  lithographj',  metallization,  and  liftoif 
steps 


Front  Metal  Lithography 

Clean  wafer  with  acetone,  methanol,  isopropyl  30  sec 

Spin  on  PMGI  at  4000  rpm  20  sec 

Bake  at  250°  C  for  5  min  on  hot  plate 

Spin  on  PMGI  at  4000  rpm  20  sec 

Bake  at  250°  C  for  5  min  on  hot  plate 

Spin  on  PMGI  at  4000  rpm  20  sec 

Bake  at  250°  C  for  5  min  on  hot  plate 

Spin  on  1400-27  lohotoresist  at  4000  rpm  30  sec 

Bake  at  100° C  for  5  min  on  hot  plate 

Edge  removal  -  expose  405  nm  UV  10  sec 

Develop  5:1  DI;351  at  500  rpm  30  sec 

DI  rinse  30  sec 

Expose  mask  (lines)  320  nm  UV  at  12  mW/cm-  for  42  sec 
Develop  5:1  DI:351  at  500  rpm  27  sec 
DI  rinse  20  sec 

DU\^  10  mW/cm-  220  nm  900  sec 
Develop  SAL  101  at  500  rpm  60  sec 
DI  rinse  30  sec 

hlask  off  structures  (if  needed)  with  1400-27  photoresist 
Bake  at  100°  C  for  5  min  on  hot  plate  if  structures  masked 

Front  Metallization 

HChDI  1:1  rinse  5  sec 
DI  rinse  30  sec 
BOE:DI  1:1  rinse  5  sec 
DI  rinse  30  sec^ 

Evaporate  200.4  titanium  sticking  layer 
Evaporate  2.5/^  gold 

Liftoff 

Acetone  soak  approximately  15  min  until  liftoff  begins 
Acetone  rinse  on  spinner  until  liftoff  complete 
Strip  PMGI  with  1165  at  low  boil  30  sec 
DI  rinse  10  sec 

Clean  wafer  with  acetone,  methanol,  isopropyl  30  sec 
Bake  at  150° C  for  6  min  on  hot  plate 

Spin  on  wax  at  2500  rpm  for  30  sec  and  put  glass  disk  on  wafer 


Figure  2.6  Expose  mask  to  photoresist-covered  GaAs  sub¬ 
strate 


resists 


Figure  2.7  Develop  first  layer  of  photoresist 

off  during  grinding  or  later  in  processing.  The  procedure  is  dangerous  to  the 
operator  because  small  GaAs  particles  become  mixed  in  with  the  slurry,  and 
GaAs  is  carcinogenic.  During  grinding  this  slurry  collects  on  various  parts  of  the 
grinding  machinery,  dries,  and  becomes  airborne  dust.  The  airborne  aluminum 
oxide  powder  is  also  hazardous.  Extreme  caution  is  required  to  avoid  inhalation 


Figure  2.9  Evaporate  titanium  and  gold  metal  layers 

of  this  material.  In  short,  this  procedure  is  fraught  with  peril. 

Table  2.6  gives  the  steps.  Figure  2.13  and  Table  2.7  show  the  measured 
thicknesses  of  the  wafers  which  sur^^ved  to  this  point. 

2.3.5  Vias  -  third  mask  level 

Yias  are  another  fabrication  step  which  is  dangerous  to  the  wafer; 


Figure  2.10  Liftoff  unwanted  metal  with  acetone 


gold 


Figure  2.11  Finished  structure  with  vertical  edge 

it  is  easA"  to  misalign  the  via  holes  and  coplanar  transition  pads,  and  it  is 
eas3‘  to  overetch.  This  step  requires  the  handling  of  hot  sulfuric  acid,  so  it  is 
dangerous  to  the  operator  as  well.  Table  2.7  lists  the  process  steps,  and  Figures 
2.14  through  2.17  show  them  in  some  detail.  There  are  no  alignment  meirks 
on  these  masks  (in  hindsight,  there  should  have  been)  because  alignment  was 
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Figure  2.12  Vertical  gold  edge  formed  with  liftoff  process 


Table  2.6  Wafer  thinning  steps 

Lapping  and  Polishing 

Complete  wafer  mounting  \\uth  vacuum  chuck 

Lap  wafer  to  125  p  with  5  ft  alimiinum  oxide/DI  water  slurry 

Polish  wafer  to  100  tolerance  with  0.3  fi 

aluminrun  oxide/DI  water/Chlorox  slurry 


accomplished  by  centering  holes  within  the  transition  pads.  The  use  of  a  mask 
aligner  with  infrared  capability  was  required  so  that  alignment  between  the  top 
and  bottom  of  the  wafer  could  be  realized. 

Figme  2.18  shows  a  pair  of  etched  via  holes  underneath  a  coplanar 
transition.  Figure  2.19  shows  a  pair  of  via  holes  from  the  bottom  side  after 
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Figure  2.13  Location.s  for  wafer  thickness  measurements 

the  thick  gold  ground  plane  has  been  electroplated  on.  Figure  2.2  earlier  in 
this  chapter  shows  the  cross  section  of  this  transition,  and  it  gives  diameters 
for  the  tops  and  bottoms  of  the  via  holes.  These  numbers  are  verj’  close  to  the 
numbers  in  the  photograph;  also  the  slightly  convex  cone  shape  of  the  vias  shown 
in  Figure  2.2  is  characteristic  of  the  chemical  etch.  If  extreme  care  and  close 
monitoring  are  exercised,  very  good  and  repeatable  via  holes  are  obtainable. 
This  step,  though,  depends  very  much  on  the  previous  thinning  step.  If  the 
wafer  is  not  of  a  verj^  uniform  thickness  (here,  “very  uniform’’  means  that  the 
wafer  thickness  h  varies  no  more  than  about  ±3  /x),  then  the  via  hole  diameters 


Table  2.7  Measured  wafer  and  ground  plane  thicknesses 
(microns) 


Wafer 

Site  1 

Site  2 

Site  3 

Site  4 

Site  5 

Ground  Plane 

1 

108 

104 

114 

105 

101 

5 

3 

115 

109 

117 

113 

111 

- 

4 

106 

107 

111 

102 

107 

7 

5 

112 

108 

116 

101 

105 

7 

6 

103 

104 

112 

105 

099 

5 

7 

097 

097 

110 

103 

099 

5 

8 

095 

095 

101 

093 

095 

7 

9 

094 

092 

103 

100 

100 

6 

10 

112 

113 

116 

103 

106 

6 

12 

103 

106 

112 

108 

104 

5 

13 

103 

104 

112 

098 

104 

4 

14 

108 

105 

115 

no 

no 

6 

15 

105 

100 

107 

103 

097 

6 

will  be  very  nonuniform  across  the  wafer.  When  the  etchant  breaks  through  the 
last  little  bit  of  GaAs  near  the  end  of  the  etch,  the  hole  diameter  increases  very 
rapidlj’.  If  the  wafer  is  nonuniform,  then  some  holes  will  have  broken  through 
and  grown  very  large  while  others  have  not  yet  broken  through.  Thus  if  the 
wafer  is  very  nonuniform,  some  holes  can  overetch  and  become  larger  than  their 
transition  pads,  while  others  have  not  yet  made  a  good  enough  breakthrough 
to  ensure  a  good  pad-to-groimd  plane  contact  when  the  back  metal  is  put  on. 

The  comments  in  Table  2.7  about  keeping  the  wafer  wet  are  very 
important.  If  the  wafer  dries  out  partially,  air  bubbles  can  form  inside  the 
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Table  2.8  Via  hole  lithography  and  etching  steps 


Via  Hole  Lithography 

Clean  wafer  with  acetone,  methanol,  isopropyl,  DI  30  sec 
Bake  at  100°C  for  6  min  on  hot  plate 
Oo  plasma,  etch  10  min 
Etch  B0E:H20  1:1  for  30  sec 
DI  rinse  30  sec 

Deposit  3000A  Si3N4  at  100°C 

Clean  wafer  with  acetone,  methanol,  isopropyl,  DI  30  sec 
Bake  at  100° C  for  6  min  on  hot  plate 

Spin  on  V^aycoat  Type  3IC  negative  resist  at  3000  rpm  for  30  sec 
Bake  at  70°  C  for  20  min  in  convection  oven 
Visually  inspect  photoresist 

Expose  mask  (via.s)  405  nm  UV  18  sec  (align  via  holes 
to  pad  centers) 

Develojo  xylene  15  sec 
Rinse  n-butyl  acetate  10  sec 
Rinse  isoiDropyl  20  sec;  blow  drj' 

Visually  inspect  via  pattern 
Oo  plasma  etch  10  min  (descum) 

Etch  B0E:H20  1:1  for  30  sec 
DI  rinse  30  sec 

Visually  inspect  etched  windows  in  Si3N4 
Strip  photoresist  with  Losolin  remover  at  low  boil 
DI  rinse 


Via  Hole  Etch 

NOTE:  etchant  solution  is  3:1:1  H2SO4  :  H2O2  :  H2O 

Pour  water  into  beaker 

Slowly  add  sulfuric  acid  to  same  beaker 

Place  beaker  into  circulating  temperature  bath:  heat  to  60°  C 

Add  peroxide  to  beaker 

Wait  15  min,  then  slowly  place  wafer  into  etchant  with  Teflon  holder 
Check  wafer  every  30  min  until  holes  are  almost  through.  Keep  the 
wafer  wet  throughout  the  process.  Do  this  bj^  lifting  the 
wafer  from  the  etchant,  rinsing  in  a  beaker  of  DI  water, 
and  putting  the  wafer  into  a  petri  dish  of  DI  water. 

When  holes  are  almost  through,  let  the  etchant  sit  in  the  bath  for 
15  min  to  reduce  its  strength 
Resume  etching.  Check  wafer  every  5  min  until  flnished 


GaAs 


Figure  2.14  Expose  via  hole  mask  to  photoresist  layer  on 
backside 

etching  via  holes.  The  surface  tension  of  these  bubbles  will  keep  the  etchant 
from  making  contact  with  the  GaAs  surface,  and  the  etch  rate  of  such  a  via  hole 
will  drop  drastically.  Either  the  hole  will  become  asymmetrical,  or  it  will  not 
be  etched  all  the  way  through.  Of  course,  during  the  etching  process  the  silicon 
nitride  is  being  etched  also.  Silicon  nitride  is  resistant  -  not  impervious  -  to 
the  sulfuric  acid  etch,  so  lengthy  periods  in  the  etchant  will  tend  to  remove  the 
silicon  nitride  mask  and  to  make  the  via  hole  bases  very  large.  Many  variables, 
including  concentration,  temperature,  and  relative  proportions  of  the  etchant 
constituents,  as  well  as  the  time  history  of  the  solution,  affect  the  etch  rate. 
Cleanliness  of  the  wafer  and  the  exposed  surface  area  also  affect  etch  rate.  The 
etching  process  must  be  monitored  very  carefully. 


resist 

Si3N4 

GaAs 

gold 

Figure  2.15  Develop  photoresist 

Si3N4 

GaAs 
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Figure  2.16  Etch  silicon  nitride  “mask”  layer  and  strip  pho¬ 
toresist 

2.3.6  Back  metallization  and  demounting 

These  are  the  final  processing  steps.  It  is  also  dangerous  to  the  wafer 
(actually,  all  processing  is  dangerous  to  GaAs  because  it  is  so  fragile!).  Care 
should  be  taken  when  electroplating  the  ground  plane  that  sharp  points  on 
the  alligator  clip  do  not  clamp  too  hard  on  the  wafer.  These  sharp  points  can 
initiate  cracks  that  propagate  freely  along  the  wafer's  cleavage  planes  when  the 
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Via 


electroplated  gold 


Figure  2.17  Etch  via  hole.s.  evaporate  gold  seed  layer,  and 
electroplate  gold  ground  plane 


Figure  2.18  Freshly  etched  ^da  holes 


wafer  is  removed  from  the  glass  disk. 

Demounting  is  perhaps  the  most  stressful  step  of  all.  If  the  wafer 
is  still  whole  at  this  point,  then  it  has  already  been  subjected  to  nontrivial 
stresses.  It  is  at  its  thinnest  and  most  fragile.  If  not  for  the  relatively  thick 


Figure  2.19  Via  holes  with  electroplated  ground  plane 


ID=9  SCM=  7801)11  VERT:  8  8 
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Figure  2.20  Dektak  plot  of  electroplated  gold  ground  plane 
thickness  of  Ifi  for  Weifer  8 
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gold  ground  plane  to  lend  it  support,  the  wafer  would  probably  not  survive 
being  removed  from  the  glass  disk.  The  -svay  this  demounting  is  performed  at 
Wright-Patterson  is  somewhat  subject  to  individual  preference  and  creativity, 
but  the  way  I  performed  it  was  to  place  the  wafer  (glass  disk  side  down)  onto  the 
hot  plate  and  to  have  a  towelette  rated  for  clean  room  use  nearby.  As  the  wafer 
heated  up  and  the  wax  softened  (about  fifteen  seconds  or  so),  I  used  a  pair  of 
forceps  to  push  the  wafer  off  the  disk  and  onto  the  towelette,  using  another  pair 
of  forceps  to  keep  the  glass  disk  stationary.  The  difficulties  with  this  process 
are  that  the  wafer  may  break  (the  object  is  to  apph'  all  force  parallel  to  the 
wafer,  laot  normal  to  it  )  and  that  the  wafer  may  adhere  to  the  towelette  if  the 
wax  is  still  soft  by  the  time  the  wafer  gets  there.  A  safer,  more  reliable  way  to 
demount  wafers  definitely  needs  to  be  developed. 

Table  2.8  lists  the  relatively  simple,  but  difficult  to  perform,  process 
steps.  Figure  2.20  shows  a  Dektak  plot  verifying  the  thickness  of  one  of  the 
wafer’s  ground  planes;  this  is  Wafer  8  with  a  thickness  of  about  7/u.  The  left 
side  of  the  plot  is  the  gold  ground  plane,  and  the  very  flat  area  on  the  right  is 
the  GaAs  wafer  surface;  this  type  of  plot  is  made  with  the  wafer  upside  down. 
This  plot  also  discloses  the  magnitude  and  character  of  the  surface  roughness 
of  electroplated  gold.  Note  that  the  roughness,  though,  is  on  the  side  of  the 
conductor  away  from  the  substrate.  Since  a  seed  layer  of  gold  was  evaporated 
onto  an  optically  flat  surface  before  electroplating,  the  substrate  to  ground 
plane  interface  is  very  smooth. 


34 


Table  2.9  Back  metallization  and  demounting  steps 


Back  Metallization 

Coat  visible  portion  of  glass  disk  with  1400-27  photoresist 
Evaporate  200A  titanium  on  back 
Evaporate  1000 A  gold  on  back 
Electroplate  5  to  7  /,/  gold  on  back 

Mounting 

Heat  wafer  on  hot  plate  at  200°C  until  wafer  detaches  from  disk 
Clean  wafer  in  petri  dishes  with  trichloroethylene,  acetone,  methanol, 
isopropyl,  and  DI  water.  Allow  to  air  dry. 

Spin  on  wax  at  2500  rpm  for  30  sec  on  silicon  wafer 

Put  silicon  wafer  on  hot  plate  at  200® C;  place  GaAs  wafer  on  it 

Clean  any  remaining  residue  from  wafer  on  spinner 


CHAPTER  III 


LOSS  MEASUREMENT  METHODOLOGY 

There  are  two  primary  methods  for  measuring  loss,  and  there  are  sev¬ 
eral  sources  of  error  associated  with  these  two  methods.  But  before  loss  mea¬ 
surements  can  be  made,  the  measurement  system  (typically  an  automatic  vector 
network  analyzer  such  as  the  Hewlett-Packard  S510B)  must  be  calibrated.  In 
the  chapter  we  explore  calibration,  loss  measurement  methods,  potential  sources 
of  error,  and  relationships  among  these  topics. 

3.1  Thru-reflect-line  (TRL)  calibration 

3.1.1  What  is  calibration? 

Calibration  is  the  technique  of  measuring  standards  whose  behaviors 
are  well  known  and  using  the  sj'stem’s  responses  to  remove  the  effects  of  sys¬ 
tematic  errors  in  the  measurement  system  [20-21]. 

Within  the  context  of  a  vector  network  analyzer  (VNA),  a  “perfect” 
network  analyzer  is  assumed  (Figure  3.1.).  The  effects  of  the  connectors  or 
transitions,  intervening  cable,  and  differences  among  the  two  ports  are  lumped 
into  the  two-port  networks  labelled  “error  boxes”  (all  calibrations  discussed 
here  will  be  understood  implicitly  to  be  two-port  calibrations).  Not  shown  in 
the  figure  are  two  other  important  corrections;  the  switch  terms  and  the  leakage 
terms.  Leakage  (or  isolation)  is  whatever  signal  gets  from  one  port  to  the  other 
when  both  ports  are  terminated  in  a  matched  load.  The  switch  (or  directivity) 
terms  arise  from  the  fact  that  within  the  “perfect  network  analyzer”  box,  there 
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are  two  reflectometers  (one  per  port).  These  reflectometers  are  switched  to 
measure  forward  and  reflected  power,  but  the  switches  are  not  matched  exactlj'. 
The  switch  terms  account  for  the  mismatch  between  these  switches.  All  these 
terms  taken  together  (four  for  the  S-parameters  in  each  error  box,  two  for  the 
switching  terms,  and  two  for  the  leakage  terms)  consitute  the  twelve-term  error 
model  used  in  the  HP  Sol  OB  VNA  (also  used  in  other  models).  The  error  box 
model  can  be  reduced  to  the  eight-term  error  model  if  the  switch  terms  are 
“ratioed  out" ,  and  isolation  is  either  omitted  or  taken  care  of  separately  [20] . 

Calibrations  differ  in  the  tyjoes  and  number  of  standards  used  to  de¬ 
duce  these  error  terms.  The  usual  method  used  for  coaxial  media  with  the  8510 
is  OSLT  (open-short-load-thru);  this  is  an  undesirable  technique  for  microstrip 
because  it  is  very  difficult  to  make  short  circuits  and  accurate  matched  loads  in 
microstrip.  The  calibration  method  used  in  this  research  is  TRL  (thru-reflect- 
line). 

3.1.2  TRL  calibration  and  standards 

TRL  calibration  [20-33]  is  an  attractive  technique  for  microstrip  me¬ 
dia  because  the  standards  required  are  a  thru  (may  be  of  nonzero  length),  a 
reflect  (may  be  a  short  or  an  open),  and  a  line  (length  depends  on  the  freciuency 
over  which  the  calibration  is  to  be  performed).  Moreover,  the  reflect  standard 
does  not  have  to  be  perfect.  This  is  the  traditional  TRL  calibration  which  is 
implemented  in  the  HP  8510B  VNA.  For  microstrip  a  thru,  an  open,  and  a  line 
are  ideal;  these  standards  are  easy  to  implement  in  this  medium.  Moreover,  an¬ 
other  advantage  to  this  method  is  that  the  standards  are  in  the  same  medium 
(that  is,  the  standards  have  the  same  characteristic  impedance  Za  and  propa¬ 
gation  constant  7  as  the  device  under  test  (DUT).  Finally,  a  great  convenience 
is  that,  since  eight  error  terms  are  needed  but  ten  are  actually  measured  (four 
for  the  thru,  four  for  the  line,  and  two  from  the  reflects),  two  more  unknowns 
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Figure  3.1  Partitioning  of  a  real  network  analyzer  system 
into  a  perfect  network  analyzer,  error  boxes,  and  the  device 
under  test  (DUX) 

can  be  determined.  These  two  tmknowns  are  taken  to  be  the  actual  reflection 
coefficient  Pl  of  the  reflect  (useful  to  determine  the  open-circuit  fringing  ca¬ 
pacitance  of  a  microstrip  open  circuit)  and  the  acutal  propagation  constant  *) 
of  the  line  (from  which  e^e  can  be  determined).  Since  the  line’s  attenuation  a 
is  the  desired  parameter  in  tliis  research,  this  approach  is  extremely  attractive. 

The  general  derivation  of  the  TRL  calibration  will  be  presented  here 
[26].  This  treatment  makes  use  of  transmission  matrices  (also  known  as  R 
matrices),  which  may  be  cascaded.  First  the  S  matrix  [5.4b]  of  the  thru  is 
measured;  it  is  converted  to  the  equivalent  J?  matrix.  If  each  connector’s  effects 
are  isolated  to  their  respective  ports,  then  this  thru  can  be  thought  of  as  two 
cascaded  R  matrices,  .4  and  B,  each  of  which  contain  their  respective  connector 
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(transition)  effects.  Note  that 

[Rab]  =  |i?.4l  ■  [Be]  (3.1) 

since  these  are  transmission  matrices. 

Next  the  line  is  connected  and  its  S  matrix  [Sale]  is  measured.  This 
is  the  same  as  [Sab]  above,  with  an  added  length  of  line  between  the  A  and  B 
matrices.  [S.4Lb]  is  converted  to  [Rale]-  noting  that 

[Rale]  =  [Ra]  •  [Rl]  •  [Re]  (3-2) 

These  two  standards  alone  provide  enough  information  to  obtain  the 
frequency-  dependent  propagation  constant  of  the  transmission  line  in  which 
calibration  is  being  performed.  To  do  so,  first  consider  the  matrix 

[X]  =  (3.3) 

T  is  fully  known  since  it  is  calculated  from  the  R  matrices,  and  hence  the 
measured  S  matrices,  of  the  line  and  thru. 

Equation  (3.3)  can  be  expanded  into 

[r]  =  [B4]  ■  [Bil .  [Bfl] .  [Bb)-'  •  IB4)-'  =  [B.,] .  |Bl]  ■  (3.4) 

The  transmission  matrix  of  the  added  length  L  of  the  line  can  be 

written: 

[Bi]=[%  (3.5) 

since  the  line  has  propagation  constant  7  =  0-!-  jfi.  Postmultiplying  (3.4)  by 
[Ra]  gives 

[T]  ■  [Ra]  =  [Ra]  •  [Rl]  (3.6) 

which  can  be  expanded  to 

Ti  T2  Rai  Ra2  _  Rai  Ra2  e~~^^  0  /„ 

Tz  r4jLi?^3  Ra4\~[Raz  Raa\'[  0  e-y^\ 
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Exi^anding  these  matrices  gives  four  equations;  eliminating  some  unknowns 


gives 


_  ToTz)  =  0; 


this  complex  equation  has  2  solutions.  Setting  e*'-^  equal  to  G\  we  can  write  the 


solution  to  (3.8)  as 


G  =  B[1±  \fD] 


The  two  quantities  represented  by  Equation  (3.9)  are 


and  IG2  they  are  found  from 


^  /Ti  +  r4\  ,  ,  /,  T^T.-To.n 

(^'T' 


(3.10) 


The  proper  solution  is  indicated  by  whichever  arg,  is  negative.  This  insures 
that  the  phase  constant  has  the  proper  sign. 


Now  the  attenuation  per  unit  length  in  dB/cm  is 

_  20logio|e-^^| 

^  £(cm) 


(3.11) 


The  effective  dielectric  constant  ere  is  also  determined  from 


dL  =  arg  = 


2-KLfy/rre 


(3.12) 


30  arg 

27rL(cm)/(GHz) 


(3.13) 


The  propagation  constant  of  the  line  is  the  desired  parameter  in  this 
research,  so  we  will  not  discuss  the  derivation  of  the  reflect  standard's  actual 
reflection  coefficient  Ti,  in  detail.  The  procedure  is  to  measure  identical  though 
unknown  and  possibly  fairly  poor  reflections  on  ports  1  and  2.  It  is  necessary 
to  know  if  the  reflect  is  an  approximate  open  or  an  approximate  short  circuit; 


this  resolves  the  ambiguitj'’  is  finding  the  correct  phase  angle  of  the  measured 
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reflection  standard’s  Fi,. 

Also,  if  so  desired,  depending  upon  the  length  of  line  in  which  the  DUX 
is  embedded  from  the  reference  plane,  be  it  positive  (extra  line)  or  negative 
(shorter  line),  the  DUT's  measured  S-parameters  may  be  de-embedded  [26]. 
This  simply  means  that  the  effects  of  the  extra  (positive  or  negative)  length  of 
line  are  removed  from  the  DUT’s  measured  S-^parameters.  and  this  is  possible 
because  we  know  both  the  physical  length  of  line  in  which  the  DUT  is  embedded, 
as  well  as  the  propagation  constant  of  said  line.  This  additional  effort  was  not 
required  in  this  research,  but  is  fairly  simple  to  accomplish. 

Now  that  we  have  the  calibration  method  in  hand,  we  need  to  examine 
standard  selection.  First,  though,  a  word  about  the  reference  plane.  The  wafer 
probe  is  an  intrusi^^e  device  which  perturbs  the  quasi-TEM  mode  fields  irear 
the  transition;  for  that  matter,  the  transition  itself  is  a  perturbation.  For  this 
reason  it  is  recommended  that  the  reference  plane  be  located  far  enough  away 
from  the  probe  tips  or  transition  (whichever  is  nearer  the  reference  plane)  that 
the  evanescent  fields  have  died  out.  Numerical  experiments  [-34]  indicate  that 
about  200  n  should  be  sufficient;  this  is  about  2  substrate  heights.  The  setup 
on  these  wafers  puts  the  reference  plane  500  jj  (oh)  away  from  the  transition; 
this  is  half  the  physical  length  of  the  non- zero-length  thru  (1000  /.i). 

Either  the  thru  or  the  reflect  can  be  used  to  determine  the  reference 
plane,  but  I  have  found  it  more  convenient  to  use  the  thru.  It  is  slightly  more 
intuitive,  and  it  is  easier  to  keep  track  of  when  probing  to  look  at  the  thru 
standard  and  note  that  the  reference  plane  is  half  of  its  length.  The  reflect 
can  extend  past  the  reference  plane,  but  if  it  extends  a  quarter-wavelength, 
then  the  open  looks  like  a  short,  so  it  is  usually  better  to  let  the  open  be  at  or 
very  slightly  past  the  reference  plane.  It  is  also  useful  to  have  a  pair  of  reflects 
facing  each  other  so  that  their  measurement  can  be  made  simultaneously  on 
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both  ports.  MMIC  fabrication  technology  is  sufficiently  advanced  that  the  two 
reflects  will  be  essentially  identical:  more  attention  should,  in  fact,  be  focussed 
on  making  identical  microwave  wafer  i^robe  contacts. 

^^’ith  the  thru  and  reflect  seen  to,  only  the  line  remains.  It  is  here 
that  confusion  is  most  likely  to  arise,  though  the  only  parameter  left  to  vary 
is  the  line  length.  TRL  requires  that  the  difference  in  line  length  between  the 
thru  and  line  be  nearly  an  odd  multiple  of  a  quarter-wavelength,  with  an  exact 
quarter-wave  being  optimal.  This  points  out  i^erhaps  the  strongest  disadvantage 
to  using  TRL:  that  is,  it  is  an  inherently  narrowband  technique.  OSLT  can  be 
made  to  cover  the  entire  100  MHz  to  40  GHz  band,  but  TRL  suffers  from  a 
bandwidth  limitation.  If  the  line  length  difference  approaches  a  multiple  of  a 
half- wavelength,  then  there  is  an  ambiguity  that  cannot  be  resolved.  It  is  known 
from  transmission  line  theory  that  adding  an  additional  half-wavelength  of  line 
to  a  transmission  line  (with  the  same  Z©,  or  course)  does  not  change  the  line 
any,  except  for  the  additional  loss.  Therefore  the  measurement  of  the  thru  and 
the  line  appear  identical  (except  for  the  aforementioned  additional  loss),  and  an 
additional  standard  will  still  be  required.  The  normal  criterion  for  selecting  the 
line  standard’s  length  is  that  the  difference  in  electrical  length  between  the  line 
and  the  thru  must  not  approach  within  20°  of  a  multiple  of  a  half-wavelength. 
That  is,  the  difference  in  phase  must  remain  between  20°  and  160°,  or  between 
200°  and  340°,  and  so  on.  Of  course,  very  long  lines  get  unwieldy  in  terms  of 
substrate  area  required  and  maintainability  of  uniformity.  The  two  equations 
used  for  calculating  the  standard’s  electrical  length  are: 


A^(deg)  = 


360°  -f-Al 


(3.14) 


and 


delay  (sec)  = 


(3.15) 
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where  /  is  the  frequency*  a,t  the  upper  or  lower  end  of  the  band.  AI  is  the 
difference  in  physical  length,  and  c  is  the  speed  of  light.  Equation  (3.14)  is 
necessary,  of  course,  to  check  that  the  electrical  length  difference  requirement  is 
not  violated,  and  this  rec^uirement  must  hold  across  the  entire  frequency  band. 
This  means  that  TRL  is  restricted  to  an  S:1  bandwidth,  unless  several  lines 
are  used  and  then  the  areas  joined  after  calibration.  This  method  is  notorious, 
however,  for  producing  jump  discontinuities  at  the  joined  band  ends,  and  is  not 
used  in  this  research.  The  6  GHz  to  40  GHz  band  is  less  than  8:1.  Equation 
(3.15)  is  required  because  the  HP  8510B  expects  the  standard  definitions  to  be 
given  in  terms  of  electrical  delay  rather  than  electrical  phase  difference.  Note 
that  specifying  a  TRL  line  standard  requires  that  a  fairly  accurate  calculation 
of  that  line’s  effective  relative  dielectric  constant  be  performed.  After  all,  there 
is  no  other  waj'  to  determine  a  line's  electrical  length  beforehand. 

Before  we  progress,  I  must  note  that  single-thru  TRL  calibration  gives 
a  fairly  noisy  measurement  of  a;  after  all,  “good”  transmission  lines  are  low-loss. 
In  addition,  I  needed  to  resolve  fairly  small  changes  in  q  due  to  the  addition 
of  grooves  of  different  orientations.  If  the  noise  in  the  extraction  of  a  were 
greater  than  the  desired  changes  in  a,  it  would  be  exceedingly  difficult  to  make 
good,  meaningful  measurements  without,  perhaps,  resorting  to  soi^histicated 
statistical  methods  based  on  the  characteristics  of  the  desired  measurements 
and  the  noise.  But  a  solution  exists. 

3.1.3  NIST  multi-line  TRL  calibration 

Dylan  Williams  and  Roger  Marks  at  the  National  Institute  of  Stan¬ 
dards  and  Technology  (NIST)  have  developed  a  very  valuable  extension  to  the 
TRL  calibration  method;  it  is  referred  to  as  multi- Hne  TRL  [35-37].  It  allows 
the  use  of  multiple  lines  to  improve  calibration  accuracy  and  bandwidth.  Us¬ 
ing  a  linearized  error  analj’^sis  of  the  TRL  calibration  method,  they  arrived  at 
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optimal,  minimum-variance  estimates  of  the  calibration  constants.  This  statis¬ 
tical  treatment  allows  great  increase  in  the  accuracy  of  the  extraction  of  the 
line’s  propagation  constant.  It  is  interesting  that,  in  their  method,  using  lossy 
lines  improves  the  calibration  accuracy.  Normallj",  as  the  difference  in  length 
between  the  thru  and  line  approaches  a  multiple  of  a  half- wavelength,  the  lin¬ 
ear  error  term  becomes  infinite  for  a  lossless  line  as  its  denominator  becomes 
infinite.  The  denominator  is  the  magnitude  of  the  sine  of  the  phase  difference. 
If  the  line  is  lossy,  the  error  term  is  less  because  the  denominator  is  the  sine  of 
a  complex  number,  and  it  is  never  zero:  thus  its  inverse  is  never  infinite.  This 
effect  is  examined  quantitatively  in  [36]. 

A  noteworthy  requirement  is  that  one  of  the  extra  lines  be  very  long, 
say  an  order  of  magnitude  longer  than  the  others.  This  greatly  increases  the 
effects  of  the  line’s  loss.  After  performing  several  on-wafer  calibration  experi¬ 
ments  with  different  combinations  of  lines  (up  to  four  short  ones  and  one  very 
long  one),  that  using  all  lines  provided  an  order  of  magnitude  less  noisy  extrac¬ 
tion  of  Q  than  the  normal,  single-line  TRL  calibration.  Leaving  out  a  single  of 
the  short  lines  degraded  the  noise  slightly;  leaving  out  two  degraded  it  much 
more  noticeably.  Using  the  four  short  lines  but  leaving  out  the  long  line  was 
almost  as  bad  as  the  single-line  method.  Therefore  I  can  confirm  their  claim 
that  more  lines  give  better  calibration,  and  that  it  is  desireable  to  have  one  very 
long  line. 

I  must  state  here  that  I  (along  with  the  rest  of  the  microwave  measure¬ 
ment  commimity)  am  greatly  indebted  to  their  extension  of  the  TRL  calibration 
technique,  and  that  I  do  not  think  I  could  have  made  these  measurements  with¬ 
out  it.  Certainly  they  would  be  much  more  noisy,  the  standard  deviations  for 
several  averaged  sets  of  measurements  would  be  much  larger,  and  the  noted 
increases  in  loss  due  to  surface  roughness  would  be  much  more  uncertain. 
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3.1.4  Sources  of  error 

The  following  are  widely  recognized  sources  of  error  in  on- wafer  mea¬ 
surements  [38]: 

1.  Degraded  Probes  and  Standards 

2.  Deviations  in  Standards 

•  Substrate  thickness 

•  Via  hole  size 

3.  Repeatability 

•  Position  (x  and  y) 

•  Pressure  (c) 

•  Angle  (rotation  of  jDrobe  in  xy  plane) 

These  items  will  now  be  addressed  individually. 

Degraded  probes  and  standards  result  after  repeated  contacts  by  the 
wafer  probes.  The  gold  on  the  transitions  gets  rubbed  off  and  even  smeared;  it  is 
a  very  ductile  metal.  The  gold  on  the  wafer  probe  tips  also  gets  worn.  This  wear 
results  in  slight  differences  among  individual  connections  to  the  same  standard 
(or  DUT,  for  that  matter). 

There  will  almost  certainly  be  deviations  in  the  standards;  even  with 
careful,  individual  attention  to  fabrication,  there  will  be  differences  from  wafer 
to  wafer  and  even  across  the  same  wafer.  The  most  noticeable  difference  is  in 
substrate  thickness,  and  this  is  most  of  the  cause  of  the  difference  in  via  hole 
size.  Due  to  the  imperfect  nature  of  the  mechanical  grinding  and  polishing 
process,  there  will  almost  always  be  a  variation  of  height  /?  across  the  two-inch- 
diameter  wafer.  Usually  this  will  be  about  ±5  //;  sometimes  it  is  a  little  better, 
and  sometimes  worse.  At  least  it  was  measured,  if  not  controlled;  little  can 
be  done  to  control  this  variation.  It  is  directly  dependent  on  the  quality  of 
grinding  that  can  be  accomplished  with  the  equipment  used.  This  variation  in 
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h  will  affect  the  diameter  of  the  \da  holes;  as  they  break  through  to  the  top 
side,  the  diameter  etch  rate  increases  dramatically.  This  is  because  the  acid  has 
ceased  to  go  through  the  Avafer  and  can  now  begin  removing  the  thin  material  at 
the  edge  of  the  newly-formed  hole.  If  the  A'ia  holes  are  drasticalh'  different  from 
standard  to  standard,  then  the  condition  of  identical  connectors  is  violated. 
Fortunately  there  can  be  a  moderate  amount  of  difference  in  via  hole  diameter 
without  changing  Aua  hole  inductance  very  much.  But  the  more  serious  effect 
is  the  change  in  the  propagation  constant.  If  h  varies,  then  all  the  lines  are  not 
identical  anymore.  Furthermore,  h  can  vary  by  b%  or  so  down  the  length  of  a 
verA’  long  line;  it  is  not  clear  AA'hat  effect  using  a  nonuniform  transmission  line 
as  a  standard  has  on  TRL  or  multiline  TRL  calibration:  this  topic  should  be 
iuA'estigated  further. 

Assuming  that  probe  tip  and  transition  erosion  and  sul^strate  thickness 
A'ariation  can  be  controlled,  then  repeatability  is  the  remaining  factor,  and  it 
is  operator  dependent.  Due  to  the  nature  of  the  chuck  and  3-axis  probe  head 
micrometer,  there  can  be  placement  errors  in  planar  (xy)  placement.  A'^ertical 
placement,  and  rotational  0  alignment.  It  is  fairly  easy  to  repeatably  place 
the  probe  tips  in  x,  y,  and  6,  but  .r  is  a  bit  more  difficult;  Appendix  B  deals 
with  this  issue  in  more  detail.  One  uses  relatiA?-e  positioning  of  the  probe  tips  on 
the  transition  for  x,  y,  and  0,  but  recj[uires  judging  the  degree  of  focus  of  the 
probe  tip  and  transition,  and  the  amount  of  lateral  movement  due  to  the  probe 
tip  contacting  the  substrate.  Errors  in  affect  the  measured  IS21I;  increased 
contact  pressure  reduces  the  apparent  loss  because  a  “better”  contact  is  made. 
It  is  possible  to  crudely  repeat  contacts  in  z  by  interactively  using  |Sii  |  and  IS22 1 
to  get  about  the  same  level  of  reflection  for  each  contact.  We  are  measuring 
lines,  which  are  supposedly  matched;  therefore  small  changes  in  the  quality  of 
contact  show  up  very  well  as  changes  in  the  small  reflection  coefficients. 
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3.2  Extraction  of  loss  from  {821 1  measurements  on  long  lines 

After  carefully  calibrating  the  measurement  system  on-wafer  as  de¬ 
scribed  in  the  last  section,  a  can  be  measured  directly.  Nothing  more  is  involved 
than  measuring  IS21I,  changing  the  sigm,  and  dividing  be  the  line  length;  this 
gives  loss  in  dB/cm  (assuming  [S21I  is  in  dB).  This  result  is  true  because,  for  a 
lossless  line,  IS21I  =  0  dB.  On  the  logarithmic  scale,  any  loss  will  show  up  as  a 
negative  transmission  coefficient.  This  number  is,  by  definition,  the  loss  of  the 
line.  All  of  the  data  labelled  as  loss  from  direct  measurement  of  lines  in  this 
dissertation  was  treated  this  way.  Also,  since  most  of  the  lines  were  measured 
several  times,  a  sample  standard  deviation  could  be  calculated  for  these  lines 
as  a  function  of  frequency. 


3.3  Extraction  of  loss  from  two-port  gap-coupled  resonator  mea¬ 
surement 

3.3.1  Extraction  of  loss  from  resonator  parameters 

The  classic  extraction  of  attenuation  constant  from  measurements  of 
two-port  gap- coupled  resonators  is  widely  [5]  used  to  measure  loss  of  microstrip 
transmission  lines.  Beginning  with  the  definition  of  resonator  Q  as  where 
U  is  the  stored  energy  and  W  is  the  power  lost  per  cycle,  and  taking  a  dis¬ 
tributed  circuit  model  of  the  transmission  line  composing  the  resonator,  it  can 
be  shown  [5]  that  ■ 

Q  =  A.  (3.16) 

2a' 

where  /3  is  the  phase  constant  of  the  line  and  a  is  the  attenuation  constant. 
Since  B  is  also  equal  to  27r/Ag,  we  see  that 


where  A 


g  -  f 


is  the  guide  wavelength.  This  equation  lumps  all  the  losses 


47 


Figure  3.2  Standard  transmission  line  model  of  a  two-port 
gap- coupled  resonator 

together  for  the  resonator,  including  the  dielectric  and  radiation  losses.  We 
stated  that  the  substrates  used  in  this  research  are  "good"  dielectrics:  just  how 
good  are  they?  We  invoke  the  often-used  expression  for  the  portion  of  the 
attenuation  constant  due  to  the  lossy  dielectric  [5]: 

Orf  =  27.3 — ^ ^  tan 6  dB /unit  length  (3. IS) 

where  tan  6  is  the  loss  tangent  for  the  substrate  material.  GaAs  has  a  loss  tan¬ 
gent  of  about  0.0006  [12].  Since  Ao  is  in  the  denominator  of  Eciuation  (3. IS). 
Orf  will  be  maximized  when  Ao  is  minimized,  or  /o  is  maximized.  So.  for  the 
dielectric  used  in  this  research,  the  portion  of  the  attenuation  constant  due  to 
the  dielectric  at  (say)  35  GHz  will  be  about  0.053  dB/cm.  Looking  ahead,  we 
see  that  at  35  GHz,  the  measured  loss  of  a  smooth  strip  is  about  O.S  dB/cm. 
So  about  6.5%  or  so  of  the  measured  loss  of  a  smooth  resonator  will  be  due  to 
the  dielectric  loss.  At  S  GHz.  the  dielectric  attenuation  is  about  0.012  dB/cm. 
whereas  the  strip  loss  is  0.37  dB/cm:  therefore  at  lower  frequencies  only  about 
3.3%  of  the  resonator  loss  is  due  to  the  dielectric.  These  percentages  are  rea¬ 
sonable  small.  The  best  approach  would  have  been  to  independently  measure 
tan  6,  but  accurately  making  this  measurement  was  not  possible. 
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Figure  3.2  is  the  standard  transmission  line  model  of  a  two-port  gap- 
coupled  resonator.  The  gaps  are  represented  here  by  equivalent  capacitive  pi- 
networks.  It  is  necessary  to  maintain  sufficient  coupling  to  excite  the  resonator 
and  to  measure  its  response,  but  light  enough  coupling  that  the  fifty-ohm  loads 
presented  by  the  network  anah'zer's  ports  do  not  unduly  perturb  the  resonator. 
It  is  possible  to  use  design  formulas  to  compute  the  shunt  and  series  capacitances 
for  this  network  and  to  calculate  losses  from  that,  but  it  is  much  easier  and  also 
sufficient  to  use  the  veiy  simple  standard  RLC  resonant  circuit  model  in  Figure 
3.3.  The  response  of  this  circuit  may  be  shown  to  be 


(3.19) 


where,  as  usual,  fo  is  the  resonant  frequency  (l/27r\/LC  in  the  circuit  model) 
and  Qo  is  the  unloaded  Q.  In  practice,  what  is  measured  is  the  loaded  Q, 
designated  Qi\  these  two  are  related  b}*  the  coupling  coefficient  k  via 


Qo  =  Q/  (1  +  ^'■) 


(3.20) 


where  the  couiDling  coefficient  is  computed  from 


k  = 


IS 


21  Imaa* 


1-|S 


21  \max 


(3.21) 


These  equations  [39]  were  used  to  calculate  the  resonator  losses  given  in  Chapter 

4. 

3.3.2  Nonlinear  curve  fit  to  standard  circuit  model 
» 

Let  us  discuss  the  fit  of  the  previously  presented  circuit  model  to  the 
measured  data.  If  there  are  too  many  second-order  effects  in  the  constructed 
resonators,  then  this  is  not  a  good  model;  we  hope  it  is  sufficient.  To  that  end, 
let  us  quickly  investigate  the  behaviors  of  the  three  parameters  that  can  be 
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Figure  3.3  Standard  series  RLC  circuit  model  for  resonator 

varied  to  fit  the  model  to  the  data. 

Figure  3.4  shows  the  effect  of  varying  IS21  |mar  on  the  model's  response; 
as  expected,  the  curve  is  shifted  vertically  Figure  3.5  shows  the  result  of  chang¬ 
ing  the  resonant  frecpiency  fo-  again,  as  expected,  the  curve  shifts  horizontally. 
Finally,  Figure  3.6  shows  the  effect  of  varying  the  parameter  Q.  This  narrows 
or  widens  the  graph,  as  one  Avould  surmise  with  /o  and  the  maximum  amplitude 
held  constant. 

Knowing  that  the  standard  RLC  circuit  model  is  a  good  model  for  res¬ 
onators.  the  ciuestion  is  how  well  it  models  the  behavior  of  the  actual  resonators 
measured  in  this  research.  Figures  3.7  and  3.8  are  actual  fitted  data  and  the 
accompanying  curves  for  what  are  quahtativelj'  termed  “poor"  and  "good"  fits, 
respectively. 

The  “poor"  fit  is  representative  of  the  worst  fit  set  of  data  taken,  and 
the  "good"  fit  is  representative  of  the  best  fit  set  of  data  taken.  As  is  usually 


50 


Figure  3.4  Effect  of  varying  IS21  \max  on  the  nonlinear  curve 
fit 


the  case  with  this  tjqje  of  measurement,  the  “average’"  fit  was  somewhat  in 
between  these  two  extremes.  Tables  of  the  fitted  parameters,  along  with  their 
residuals,  are  presented  in  Chapter  4.  One  observation:  for  the  “poor”  fits,  one 
side  of  the  model  was  below  the  data,  while  the  other  side  was  above  the  data. 
This  behavior  was  fairly  consistent  when  the  fit  was  “poor,”  but  the  cause  is 
unknown  at  this  time.  Also  the  term  “poor”  fit  is  relative.  Very  good  agreement 
between  the  fitted  model  and  the  actual  data,  was  observed  in  everj^  case  but 
two  out  of  84.  When  the  standard  deviations  were  calculated,  these  two  q’s 
were  more  than  one  standard  deviation  from  the  mean,  so  they  were  discarded 
and  the  means  were  calculated' again  without  them,  as  is  standard  statistical 
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Figure  3.5  Effect  of  varj'ing  /©  on  the  nonlinear  curve  fit 

practice. 

Other  investigators  have  examined  fitting  this  kind  of  model  to  their 
data.  [40].  The  model  was  fit  in  this  case  by  minimizing  the  sum  of  the  squared 
errors  given  by  the  following  equation: 


-1^21  \maT 


(3.21) 


By  the  way.  this  is  obviously  a  nonlinear  curve  fit.  To  do  a  standard  least- 
squares  curve  fit  would  involve  finding  the  values  of  the  three  parameters  for 
which  the  following  are  true: 


d\S2i\  max 


=  0,  and 


(3.22) 


Figure  3.6  Effect  of  varying  Qi  on  the  nonlinear  curve  fit 

As  the  conditions  in  (3.22)  give  rise  to  some  very  nasty-looking  equations  (which 
are  not  reproduced  here  to  show  mercy  to  the  reader),  these  fits  were  done 
pseudo-manually.  That  means  that  the  data  were  imported  into  Borland's 
Quattro^-^^  Pro  4.0  and  the  nonlinear  optimizer  therein  was  used  to  minimize 
Equation  (3.21).  But  the  optimizer  had  a  hard  time  finding  the  minimum; 
in  fact,  it  was  reluctant  to  vary  Qi  much  at  all.  Therefore  manual  tweaking- 
varying  a  parameter  and  noting  its  effect  on  the  sum  of  squared  errors-was 
used,  first  on  the  resonant  frequency  and  maximum  amplitude,  and  then  on 
the  Qi,  with  some  going  back  and  forth.  Generally  it  was  easy  to  find  /o,  and 
the  optimizer  had  no  problem  with  that.  Maximum  amplitude  was  somewhat 
harder  to  find,  but  the  optimizer  seemed  to  handle  that  all  right,  also.  The 
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■  Data  - Model 

Figure  3.7  Example  of  a  ‘"poor"  fit 

loaded  quality  factor  required  frequent  human  intervention.  It  should  be  said 
that  all  this  fitting  was  necessary  because  the  data  were  verj?’  sparse.  If  the 
data  were  very  dense,  then  the  resonant  frequency  and  maximum  amplitude 
could  have  been  read  right  off  the  graph,  along  with  the  -3  dB  points,  and 
Qi  could  have  been  calculated  the  old-fashioned  way,  as  resonant  frequency 
divided  by  bandwidth.  However,  a  broadband  calibration  was  performed,  and 
the  resonators  were  used  in  several  modes;  thus  data  were  sparse.  Only  about 
five  data  points  per  resonance  were  above  the  -3  dB  points;  for  each  peak 
eleven  data  points  were  used  in  the  final  parameter  extractions.  Beyond  this 
point  noise  dominates  the  data,  and  the  model  should  only  be  very  accurate  in 
the  vicinity  of  the  resonance  anyway. 


IS21I 
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Figure  3.8  Example  of  a  ‘'good’*  fit 

3.4  Consideration  of  ground  plane  loss 

Holloway  [4]  gives  in  his  PhD  dissertation  an  expression  for  the  atten¬ 
uation  constant  due  to  the  loss  in  the  ground  plane  of  a  inicrostrip  of  width  iv 
and  height  h: 


is  the  usual  surface  resistance  in  ohms/square.  fi  and  a  are  the  permeability  and 
conductivity  of  the  ground  plane  metal,  respectively.  This  expression  assumes  a 


55 


- 96u  . 98u  - 100u 

Figure  3.9  Numerically  computed  ground  plane  loss 

constant  strip  current  distribution  (current  density  equals  the  total  strip  current 
divided  by  the  strip  width).  Holloway  used  a  standard  wall  loss  perturbation 
analysis;  he  also  neglected  the  transverse  component  of  the  current  and  worked 
only  with  .r-directed  magnetic  vector  potential. 

The  expression  in  (3.23)  is  integrated  numerically  to  obtain  ground 
plane  losses  which  are  extracted  from  total  losses  (in  direct  S21  measurements) 
in  order  to  give  the  total  loss  of  the  strip  and  that  added  by  the  surface  rough¬ 
nesses.  Figure  3.9  shows  this  ground  plane  loss  for  a  gold  ground  plane  with 
73|U-wide  strips  on  GaAs  substrates  of  varying  height  h. 


CHAPTER  IV 


LOSS  MEASUREMENT  RESULTS  AND  COMPARISONS 

The  measurement  results  are  summarized  and  presented  in  this  chap¬ 
ter.  Not  all  fifteen  wafers  survived,  of  course,  so  less  data  were  gathered  than 
was  desired.  Nevertheless,  a  large  amount  of  information  was  amassed,  and 
enough  was  taken  to  perform  some  statistical  analysis  and  to  make  some  con¬ 
clusions.  It  should  be  noted  at  this  point  that  whenever  reference  is  made  in 
this  chapter  to  the  mean  or  standard  deviation,  that  it  is  the  sample  mean 
and  sample  standard  deviation  which  are  being  mentioned.  The  practical  con¬ 
sequence  of  this  distinction  is  that  the  denominator  of  the  calculation  of  the 
standard  variance  5^  is  iV  —  1  instead  of  N.  Since  the  mean  is  calculated  from 
a  limited  number  of  observations  and  is  not  known  a  'priori,  using  N  would 
underestimate  the  standard  deviation  s  [41-42]. 

4.1  TRL  extraction  results 

We  have  seen  that  effective  dielectric  constant  and  loss  can  be  ex¬ 
tracted  from  the  TRL  calibration  process.  The  only  such  data  presented  here, 
though,  are  for  the  smooth  lines.  No  standards  were  made  with  parallel  or 
transverse  grooves  (in  retrospect,  a  regrettable  omission).  Of  course,  with  zero- 
length  thru  and  reflect  standards,  only  the  lines  would  need  grooves.  The  long 
lines  obviously  are  grooved,  but  none  of  the  short  lines  are.  This  smooth  line 
data,  though,  can  be  compared  with  IS21 1  measurements  and  resonator  data  for 
smooth  lines. 

Figure  4.1  is  a  graph  of  the  effective  relative  dielectric  constant  e^e 
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Figure  4.1  Effective  relative  dielectric  constant  extracted 
from  TRL  calibration  on  wafer  6 

of  the  lines  on  wafer  6.  The  solid  line  is  an  average  of  the  ten  measurements; 
the  dashed  lines  are  e^e  ±5.  Note  that  ere  increases  with  frequency  as  theory 
predicts:  as  frequency  increases,  more  of  the  electric  field  lines  are  confined 
underneath  the  strip  conductor  and  thus  in  the  dielectric.  Therefore  tre  rises. 
Note  also  that  the  standard  deviation  s  of  these  measurements  is  roughlj"  0.004 
across  the  band,  indicating  a  high  degree  of  certainty  in  the  measured  values  of 

€re . 

Figure  4.2  is  a  graph  of  the  loss  in  dB/cm  from  the  same  ten  measure¬ 
ments  on  wafer  6.  Like  Figure  4.1,  the  solid  line  is  the  mean  and  the  dashed 
lines  are  one  standard  deviation  on  each  side  of  the  mean.  For  graphs  of  single 
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Frequency  (GHz) 


Figure  4.2  Loss  extracted  from  TRL  calibration  on  wafer 
6 

parameters  with  one  solid  and  two  dashed  lines  in  the  rest  of  this  dissertation, 
the  same  meaning  will  apply  (solid  hue  is  computed  or  measured  value  and 
dashed  lines  are  value  plus  or  minus  one  standard  deviation).  As  expected  the 
loss  increases  with  frequency. 

4.2  Direct  |S2i|  results  from  long  lines 

Here  we  examine  comparisons  between  lines  with  different  surface 
roughnesses.  The  comparisons  are  somewhat  misleading  at  first,  because  these 
next  several  graphs  have  not  had  the  ground  plane  loss  extracted  yet.  Therefore 
thej^  show  the  effects  of  a  rough  surface,  but  only  the  strip  is  rough.  Of  course, 
the  strip  dimensions  were  chosen  so  that  most  of  the  loss  would  be  contributed 


59 


Figure  4.3  Smooth  loss  extracted  from  line  2  [Soi  |  on  wafer 
6 

by  it,  anyway. 

Figures  4.3,  4.4,  and  4.5  show  the  total  loss  due  to  lines  with  no 
grooves,  parallel  grooves,  and  transverse  grooves,  respectively.  As  before,  these 
are  means  of  ten  sets  of  measurements,  with  standard  deviations  indicated  by 
dotted  lines.  Figure  4.6  shows  the  means  all  on  the  same  graph  so  that  they 
may  be  compared.  Figure  4.7  is  the  same  graph  at  Figure  4.6  but  is  strip  loss 
only  (ground  plane  losses  have  been  subtracted). 

4.3  Two-port  gap-coupled  resonator  results 

Resonator  data  exist  from  wafers  5  and  6.  Smooth,  parallel,  and  trans¬ 
verse  resonator  data  exist  from  wafer  5,  and  these  data  are  presented  in  Figure 
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Figure  4.4  Parallel  loss  extracted  from  line  3  |S9i  |  on  wafer 
6 

4.8.  Only  four  frequencies  are  given  because  only  one  (length)  resonator  is 
on  this  wafer,  and  it  has  four  resonances  in  the  six-to-forty  gigahertz  band. 
The  lower  two  resonances  give  losses  that  fall  on  top  of  each  other  for  smooth 
and  parallel  losses.  The  losses  in  Figure  4.9  are  from  wafer  6,  but  no  smooth 
resonator  survived;  also  different  resonator  lengths  are  on  board.  This  means 
that  some  loss  points  are  not  directly  comparable  because  they  are  at  different 
frequencies. 

All  of  the  resonator  nonlinear  curve  fit  statistics  are  in  Tables  4.1 
through  4.3.  Only  one  set  of  measurements  of  RIP  was  taken;  this  particular 
structure  may  be  partially  damaged.  Only  two  sets  of  measurements  of  R3S, 
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Frequency  (GHz) 


Figure  4.5  Transverse  loss  extracted  from  line  2  IS21I  on 
wafer  6 

R3P,  and  R3T  from  Wafer  5  were  taken  due  to  the  need  of  the  probe  station 
for  other  measurements. 

4.4  Data  from  other  methods  and  sources 

Other  data  are  now  presented.  The  following  numbers  are  those  ar¬ 
rived  at  by  other,  more  (but  sometimes  less)  traditional  methods.  The  more 
traditional  methods  include  approximate  calculations  by  design  equations  em¬ 
pirically  arrived  at,  and  numerical  results  from  commercial  simulation  software. 

4.4.1  Approximate  calculations 

One  of  the  latest  refinements  of  the  standard  approximate  design  for- 
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Frequency  (GHz) 


Figure  4.6  Comparison  of  smooth,  parallel,  and  transverse 
total  losses 


mulas  is  presented  below.  These  equations  are  garnered  from  other  sources  and 
presented  verbatim  so  that  they  do  not  have  to  be  looked  up,  ans  to  show  what 
equations  were  actually  used.  These  equations  (or  similar  forms)  also  seem  to 
be  undergoing  constant  revision  to  improve  their  accuracies  over  desired  ranges 
of  parameters. 

Two  numbers  made  use  of  in  these  formulas  are  the  surface  resistivity 


R$ 


(f] /square) 


(4.1) 


where  a  is  the  conductivity  of  the  metal  and  6  is  the  familiar  skin  depth 
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Frequency  (GHz) 


Figure  4.7  Comparison  of  smooth,  parallel,  and  transverse 
losses  of  strip 


where  uj  is  the  radian  frequenc}'  of  operation  and  ^  is  the  permeability  of  the 
metal.  In  this  dissertation  fx  is  taken  to  be  equal  to  Ho  since  onlj'  nonferrous 
metals  (i.e.,  gold)  are  considered. 

The  effective  dielectric  constant  Cre  is  given  by 


^  er  -  1 
9  O 


and 


<  1),  (4.4) 


where  Cr,  w,  and  h  are  as  previously  defined  (substrate  relative  dielectric  con- 


Loss  (dB/cm) 


64 


Figure  4.8  Comparison  of  smooth,  parallel,  and  transverse 
resonator  losses  on  Wafer  5 


stant,  strip  width,  and  substrate  height,  respectively). 

The  characteristic  impedance  Zo  is  found  from 


in  which  jyo  =  is  the  characteristic  impedance  of  free  space. 

All  these  are  quasistatic  (/  =  0)  quantities  for  zero-thickness  {i  =  0) 


microstrip.  There  exist  corrections  for  nonzero-thickness  strips  and  for  dis- 
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Frequency  (GHz) 


*  Parallel  x  Transverse 


Figure  4.9  Comparison  of  parallel  and  transverse  resonator 
losses  on  Wafer  6 


persion  (frequency  dependence  of  e^e).  The  usual  nonzero-thickness  correction 
approach  is  to  modify  the  value  of  the  strip  width  lo  by  adding  an  additional 
incremental  width  Aw  to  arrive  at  a  corrected  width  lo'  as  follows: 


w'  =  w  +  Aw' 


(4.8) 
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Table  4.1  Wafer  6  R5P  and  RIP  Curve  Fit  Statistics 


Peak  Number 

Frequency  (GHz) 

Qo 

Residual 

a  (dB/cm) 

R5P 

1 

12.8184 

67.7785 

3.11E-6 

0.4980 

1 

12.8171 

68.2594 

3.97E-6 

0.4947 

1 

12.8168 

67.9720 

3.69E-6 

0.4968 

1 

12.8173 

67.9426 

2.88E-6 

0.4970 

1 

12.8193 

67.9534 

3.65E-6 

0.4967 

1 

12.8177 

67.9577 

3.74E-6 

0.4969 

1 

12.8177 

67.9740 

4.44E-6 

0.4967 

1 

12.8154 

67.9490 

1.63E-6 

0.4967 

a,vg 

12.8175 

- 

- 

0.4967 

std 

00.0011 

- 

- 

0.0009 

2 

25.6001 

84.2145 

1.07E-6 

0.8030 

2 

25.5969 

83.6461 

3.82E-6 

0.8088 

0 

25.5967 

84.6621 

2.78E-6 

0.7991 

2 

25.5987 

84.7382 

2.68E-6 

0.7985 

2 

25.6015 

84.5925 

1.59E-6 

0.7995 

2 

25.5984 

84.5951 

2.49E-6 

0.7999 

2 

25.5980 

84.6033 

3.37E-6 

0.7997 

2 

25.5951 

84.4161 

7.27E-6 

0.8013 

avg 

25.5984 

- 

- 

0.8001 

std 

00.0019 

- 

- 

0.0014 

3 

38.3189 

87.3254 

5.53E-6 

1.1647 

3 

38.3135 

89.4877 

1.55E-5 

1.1370 

3 

38.3082 

87.1952 

5.12E-6 

1.1668 

3 

38.3101 

87.2507 

9.26E-6 

1.1661 

3 

38.3141 

87.3270 

6.13E-6 

1.1646 

3 

38.3095 

87.3416 

1.40E-5 

1.1650 

3 

38.3093 

87.3979 

1.50E-5 

1.1640 

3 

38.3073 

87.1123 

6.81E-6 

1.1678 

avg 

38.3113 

- 

- 

1.1657 

std 

00.0040 

- 

- 

0.0013 

RIP 

1 

06.9232 

52.8352 

1.29E-6 

0.3450 

9 

13.8371 

72.2057 

7.49E-6 

0.5047 

3 

20.7052 

87.6662 

2.12E-6 

0.6229 

4 

27.5638 

90.1732 

3.94E-5 

0.8079 

5 

34.3611 

98.6639 

5.00E-5 

0.9226 
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Table  4.2  Wafer  6  R2T  Curve  Fit  Statistics 


Peak  Number 

Frequency  ( GHz) 

Qo 

Residual 

a  (dB/cm) 

1 

07.9054 

52.0458 

1.92E-6 

0.3998 

1 

07.9054 

52.0390 

2.02E-6 

0.4001 

1 

07.9033 

51.0160 

2.10E-6 

0.4080 

1 

07.9033 

51.0228 

1.53E-6 

0.4080 

1 

07.9033 

51.0186 

1.90E-6 

0.4078 

1 

07.9034 

51.0118 

2.19E-6 

0.4082 

1 

07.9034 

51.0179 

2.32E-6 

0.4080 

1 

07.9034 

51.0000 

1.84E-6 

0.4083 

avg 

07.9034 

- 

- 

0.4081 

std 

5E-5 

- 

- 

0.0002 

2 

15.7946 

75.1066 

l.llE-5 

0.5540 

2 

15.7921 

75.0424 

1.22E-5 

0.5546 

2 

15.7913 

75.1256 

1.28E-5 

0.5540 

2 

15.7930 

75.0975 

1.19E-5 

0.5543 

2 

15.7948 

75.1144 

1.33E-5 

0.5539 

2 

15.7915 

75.1134 

1.23E-5 

0.5542 

2 

15.7921 

75.1118 

1.24E-5 

0.5541 

2 

15.7921 

75.1550 

1.14E-5 

0.5539 

a.vg 

15.7927 

- 

- 

0.5541 

std 

00.0013 

- 

- 

0.0002 

3 

23.6545 

80.9783 

3.58E-5 

0.7711 

3 

23.6448 

84.1463 

7.22E-5 

0.7421 

3 

23.6473 

85.4495 

8.13E-5 

0.7309 

3 

23.6527 

86.4273 

9.02E-5 

0.7228 

3 

23.6489 

85.3974 

2.24E-5 

0.7310 

3 

23.6463 

85.3704 

2.79E-5 

0.7317 

3 

23.6476 

85.4406 

2.22E-5 

0.7310 

3 

23.6310 

85.4213 

5.20E-5 

0.7304 

avg 

23.6442 

- 

- 

0.7310 

std 

00.0067 

- 

- 

0.0004 

4 

31.4163 

88.2848 

1.08E-5 

0.9418 

4 

31.4114 

88.2899 

1.90E-5 

0.9421 

4 

31.4098 

88.1432 

1.02E-5 

0.9437 

4 

31.4105 

88.1154 

8.15E-6 

0.9440 

4 

31.4157 

88.1831 

8.19E-6 

0.9429 

4 

31.4099 

88.1779 

1.21E-5 

0.9433 

4 

31.4099 

88.1278 

8.47E-6 

0.9437 

4 

31.4099 

76.7000 

9.60E-6 

0.9447 

avg 

31.4110 

- 

- 

0.9435 

std 

00.0020 

- 

.  - 

0.0008 

5 

39.1577 

93.3466 

2.26E-5 

1.1139 

5 

39.1497 

91.3018 

1.19E-5 

1.1391 

5 

39.1476 

93.4342 

2.48E-5 

1.1131 

5 

39.1510 

93.5007 

1.97E-5 

1.1124 

5 

39.1554 

93.5223 

2.22E-5 

1.1117 

5 

39.1502 

93.4578 

2.54E-5 

1.1130 

5 

39.1484 

93.4946 

7.08E-5 

1.1123 

5 

39.1473 

93.6408 

1.38E-5 

1.1106 

avg 

39.1511 

- 

- 

1.1124 

std 

00.0037 

- 

- 

0.0010 

Table  4.3  Wafer  5  R3S,  R3P,  R3T  (Resonator  3:  Smooth, 
Parallel,  Transverse)  Curve  Fit  Statistics 
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Peak  Number 

Frequency  (GHz) 

Qo 

Residual 

Q  (dB/cm) 

1 

09.3879 

R3S 

67.4902 

9.86E-6 

0.3651 

1 

09.3898 

63.2882 

7.46E-6 

0.3891 

avg 

09.3889 

- 

- 

0.3771 

2 

18.7351 

81.2699 

1.56E-5 

0.6062 

2 

18.7369 

81.2633 

1.64E-5 

0.6060 

avg 

18.7360 

- 

- 

0.6061 

3 

28.0024 

94.8330 

4.01E-6 

0.7787 

3 

28.0053 

93.6806 

1.82E-5 

0.7879 

avg 

28.0039 

- 

- 

0.7833 

4 

37.1780 

93.8788 

6.53E-6 

1.0489 

4 

37.1841 

89.4470 

1.65E-5 

1.1005 

avg 

37.1811 

- 

- 

1.0747 

1 

09.3573 

R3T 

68.7296 

1.36E-5 

0.3573 

1 

09.3587 

63.4497 

9.19E-6 

0.3869 

avg 

09.3580 

- 

- 

0.3721 

2 

18.6921 

81.0044 

9.04E-6 

0.6068 

2 

18.6951 

80.9543 

6.06E-6 

0.6069 

avg 

18.6936 

- 

- 

0.6069 

3 

27.9491 

92.5263 

4.18E-6 

0.7965 

3 

27.9550 

92.4895 

5.80E-6 

0.7966 

avg 

27.9521 

- 

- 

0.7966 

4 

37.1230 

90.9298 

4.62E-6 

1.0813 

4 

37.1303 

86.6546 

5.69E-6 

1.1343 

avg 

37.1267 

- 

- 

1.1078 

1 

09.3346 

R3T 

62.0335 

9.74E-6 

0.3949 

1 

09.3339 

55.7489 

7.55E-6 

0.4391 

avg 

09.3343 

- 

- 

0.4170 

2 

18.6208 

76.5255 

2.34E-5 

0.6398 

2 

18.6230 

74.2342 

1.18E-5 

0.6593 

avg 

18.6219 

- 

- 

0.6496 

3 

27.8535 

87.0046 

2.30E-5 

0.8442 

3 

27.8587 

86.9485 

1.73E-5 

0.8444 

avg 

27.8561 

- 

- 

0.8443 

4 

36.9835 

82.4678 

1.74E-5 

1.1878 

4 

36.9800 

83.6923 

1.03E-2 

1.1696 

avg 

36.9818 

- 

- 

1.1787 
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This  U’’  is  then  substituted  back  into  the  previous  equations  in  place  of  ■«’  to 
arrive  at  the  new  values  of  Cre  and  Zo  for  t  0. 

The  dispersion  correction  is  incorporated  as  follows,  with  /  in  GHz 


and  h  in  cm: 


^re  if)  — 


6r  ^reif  —  0) 

l  +  Pif) 


P{f)  =  PiPoimSU  +  P3P4)  (10/  h)] 


(4.10) 


0  **!!?/'  /  ' 

Pi  =  0.274SS  +  0.6.315  +  - - ,  ,,,  -  -  0.0656S3e-®-"°^®“’/*  (4.11) 

(1  +  O.lo I  j  j  h 


Pi  =  0.33622  [1  -  e-°’ 


-0.03442  «r 


(4.12) 


Ps  =  0.0363  [l  -  ''  (4.13) 


P4  =  1  +  2.751  1 


(4.14) 


The  accuracj'  of  this  dispersion  correction  is  given  as  better  than  0.6%  for 


0.1  <  wjh  <  100 


1  <  <  20 


0<  /i/A<  0.13 


Effective  Relative  Dielectric  Constant 
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20  25 

Frequency  (GHz) 


h=96u 


h=98u 


h=100u 


Figure  4.10  ^reif)  calculated  from  design  formulas 


One  set  of  equations  given  for  calculating  conductor  losses  due  to  the 
skin  effect  and  finite  conductivity  is  cited  by  Wadell  [3].  For  w/h  <  1, 


TT  InlO 


hZo  e^o/60 


For  IV /h  >  1, 


Zo  Rs 
20  TT^  /i  In  10 


w  d  w 

^  ('  +  1  + IT 


(4.16) 


These  values  of  etc  are  in  dB/unit  length,  and  ^  is  defined  as  follows: 


dw  f  l  \  .  4  izw 

ar  =  Uj*"-r- 


for  1 

h  2  TT 


(4.17) 


in  which  ac,o  is  the  smooth  loss  calculated  from  Equations  (4.16)  through  (4.19), 
A  is  the  rms  surface  roughness,  and  S  is  the  skin  depth. 

These  equations  were  used  to  calculate  microstrip  parameters  for  the 
real  lines  on  the  wafers  fabricated  for  this  research.  Figure  4.10  is  the  disper¬ 
sive  effective  relative  dielectric  constant.  Figures  4.11  through  4.14  show  the 
different  variations  of  the  line  loss.  Figure  4.11  plots  a  for  h  =  96/u  with  rms 
surface  roughnesses  of  0//,  0.125ju,  and  0.25/u  (no  grooves,  0.25//-deep  grooves, 
and  0.5/y-deep  grooves,  respectively).  Figure  4.12  is  the  same  for  h  =  98//,  and 
Figure  4.13  lets  h  =  100//.  These  different  data  are  presented  because,  as  noted 
in  Chapter  2,  h  varies  across  the  wafer.  Figure  4.14  shows  how  a  varies  with 
respect  to  h  for  an  rms  surface  roughness  of  0.125//  (the  most  likely  value  of 
surface  roughness  for  Wafer  6,  based  on  scanning  electron  microscope  measure¬ 
ments).  Note  that  Figure  4.14  only  shows  the  35-40  GHz  range;  this  is  because 
the  approximate  formula  given  actually  has  a  very  small  roughness  correction 
factor.  Only  a  small  frequency  range  was  presented  so  that  the  change  in  a 
with  respect  to  rms  surface  roughness  could  be  seen. 
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- rough=Ou  - rough=0.125u - rough=0.25u 

Figure  4.11  Loss  calculated  from  design  formulas  for  h  — 

96//,  variable  roughness 

4.4.2  Hewlett-Packard  Microwave  Design  System  (MDS) 

The  commercially  available  program  Microwave  Design  System  (MDS) 
by  Hewlett-Packard  accounts  for  microstrip  loss  due  to  lossy  dielectric  (not  con¬ 
sidered  in  this  work),  finite  conductivity,  and  surface  roughness.  MDS  essen¬ 
tially  incorporates  the  phenomenological  loss  model  as  described  by  Itoh  [50]. 
The  surface  roughness  is  implemented  empirically  from  experimental  data  in  a 
very  simple  fashion. 

MDS  was  used  to  model  the  lines  as  fabricated  on  the  real  wafers. 
Nonuniform  wafer  thickness  is  not  supported  in  MDS,  so  the  simulations  were 
run  for  several  values  of  h  and  rms  surface  roughness  as  for  the  approximate 


- rough=Ou  . rough=0.125u - rough=0.25u 

Figure  4.12  Loss  calculated  from  design  formulas  for  h  = 

98//,  variable  roughness 

formulas;  this  makes  comparisons  relatively  easy.  Results  for  h  =  98//  and  a  = 
0.125//  are  presented,  as  for  the  approximate  formulas,  since  those  dimensions 
most  closely  match  what  was  actually  fabricated. 

Figure  4.15  shows  a  for  h  =  96//  with  rms  surface  roughnesses  of 
0//,  0.125//,  and  0.25//  (no  grooves,  0.25//-deep  grooves,  and  0.5//-deep  grooves, 
respectively).  Figure  4.16  is  the  same  for  h  =  98//,  and  Figure  4.17  has  h  = 
100//.  noted  in  Chapter  2,  h  varies  across  the  wafer.  Figure  4.14  shows  how 
a  varies  with  respect  to  h  for  an  rms  surface  roughness  of  0.125//  (again,  the 
most  likely  value  of  surface  roughness  for  Wafer  6)  Note  that,  like  Figure  4.14, 
Figure  4.18  only  shows  the  35-40  GHz  range  due  to  the  small  variation  in  the 


Figure  4.13  Loss  calculated  from  design  formulas  for  h  — 

100/i,  variable  roughness 

roughness  correction  at  this  small  value  of  surface  roughness. 

4.4.3  Other  experimental  data 

One  set  of  experimental  data  available  (and  often  quoted)  in  the  lit¬ 
erature  is  Goldfarb’s  and  Platzker’s  [30]  resonator  data.  They  made  careful 
loss  measurements  on  microstrip  resonators  of  varying  widths  on  a  nominal 
100/u  thick  GaAs  substrate.  These  structures  have  no  roughness;  at  least  none 
is  mentioned.  The  nearest  similar  structure  in  their  work  is  a  70yu  wide  line. 
All  the  lines  in  this  research  are  73//  wide,  so  the  results  are  similar,  but  not 
directly  comparable.  In  fact,  their  lines  should  be  slightly  more  lossy  since 
they  are  narrower.  These  data  have  been  read  off  of  a  graph  in  Goldfarb’s  and 
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- h=96u  . h=98u  - h=100u 

Figure  4.14  Loss  calculated  from  design  formulas  for  rms 
roughness  =  0.125yu,  variable  h 

Platzker’s  paper,  so  the  numbers  may  differ  slightly  from  their  actual  values; 
however,  they  are  close  enough  for  rough  comparisons.  The  data  are  presented 
in  tabular  form  here  and  are  included  on  some  of  the  graphs  presented  later. 
Note  that  the  error  bars  are  not  plotted,  but  are  tabulated. 

4.5  Comparisons 

Finally,  all  the  losses  are  brought  together.  First  of  all,  though,  Figure 
4.19  is  a  comparison  of  a  extracted  from  the  TRL  measurement  and  a  from 
S-parameter  measurements;  as  before,  these  are  actually  the  means  of  ten  sets 
of  measurements. 

The  next  three  figures  capture,  essentially,  the  essence  of  Chapters  2,3, 
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- rough=Ou  . rough=0.125u - rough=0.25u 

Figure  4.15  Loss  calculated  with  MDS  for  h  =  96^,  variable 
roughness 

and  4,  and  indeed  all  of  the  experimental  work.  Figures  4.20  through  4.22  bring 
together  the  losses  as  measured  and  calculated  in  different  ways  together  for 
comparison.  Figure  4.20  is  the  data  for  smooth  (unroughened)  microstrip.  On 
this  graph  are  S-parameter  measurements,  approximate  formula  calculations, 
MDS  calculations,  smooth  resonator  data  from  a  Wafer  5  (the  line  data  are  from 
Wafer  6),  and  Goldfarb’s  resonator  data.  As  mentioned  before,  Goldfarb's  data 
are  for  slightly  narrower  lines. 

Figure  4.21  is  contains  the  same  information,  except  that  Goldfarb 
did  not  make  measurements  on  roughened  surfaces;  therefore  his  data  are  only 
on  Figure  4.20.  Instead  resonator  data  from  Wafer  6  are  included.  Note  that 


- rough=Ou  . rough=0.125u - rough=0.25u 

Figure  4.16  Loss  calculated  with  MDS  for  h  =  98/^,  variable 
roughness 

MDS  and  the  approximate  formulas  do  not  distinguish  between  parallel  and 
transverse  grooves,  so  these  data  are  the  same  as  those  presented  in  Figure 
4.22.  That  figure  contains  the  same  information,  except  the  resonator  data  and 
S-parameter  data  are  from  structures  with  transverse  grooves.  Only  this  work 
provides  the  measurements  on  roughnesses  of  different  “polarization.” 

A  few  generalities  can  be  seen.  First,  the  approximate  formulas  seem 
to  form  an  upper  bound,  while  MDS  seems  to  give  a  lower  bound  on  the  loss. 
Line  measurements  coincide  with  MDS  or  lie  between  the  extremes,  whereas 
most  of  the  resonator  data  seems  to  be  scattered  within  these  two  bounds. 

More  specifically,  in  Figure  4.20  Goldfarb’s  data  lie  between  the  MDS 


- rough=Ou  - rough=0,125u - rough=0.25u 

Figure  4.17  Loss  calculated  with  MDS  for  h  =  100/i,  vari¬ 
able  roughness 

and  measured  IS21I  data,  corresponding  sometimes  to  one,  sometimes  to  the 
other,  and  sometimes  in  between.  The  R3S  resonator  data  agree  more  or  less 
with  the  IS21I  data,  but  the  high-frequencj'  point  is  much  higher;  in  fact,  it 
''  agrees  with  the  approximate  calculations.  Perhaps  the  resonator  loss  at  that 
point  receives  a  higher  contribution  from  dielectric  loss. 

Figure  4.21,  comparing  the  parallel  data  along  with  MDS  and  the  ap¬ 
proximate  formulas,  is  even  more  interesting.  The  measured  data  agrees  very 
well  with  the  MDS  numbers.  The  parallel  resonator  data  agree  very  well  with 
the  IS21I  data,  again  except  at  the  higher  frequencies.  R.3P  measurements, 
though  from  a  different  wafer  than  the  line  measurements,  seem  to  agree  par- 
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Frequency  (GHz) 

- h=96u  . h=98u  - h=100u 

Figure  4.18  Loss  calculated  with  MDS  for  rms  roughness 

0.125/x,  variable  h 

ticularly  well.  RoP  data,  from  the  same  wafer  as  the  line  measurements,  agree 
almost  as  well.  Not  just  one  set  of  resonator  points,  but  both  sets,  diverge  from 
MDS  and  line  measurements  at  the  high  end  and  approach  the  approximate 
formulas. 

Finally,  several  items  emerge  after  scrutinizing  Figure  4.22,  especially 
keeping  in  mind  the  previous  figure.  Except  for  the  conspicuous  dip  in  the  21-25 
GHz  or  so  range,  |S2ii  measurements  are  slightly  above  the  MDS  predictions; 
this  is  the  same  as  saying  that  they  are  generally  greater  than  the  parallel 
groove  measurements.  This  is  expected.  What  is  not  expected  is  the  dip.  Once 
more,  resonator  measurements  display  similar  behavior  to  those  in  Figure  4.21; 


so 


Table  4.4  Goldfarb  and  Platzker  Resonator  Data 


Frequency  (GHz) 

a  (dB/cm) 

Error  Bar  (dB/cm) 

07.5 

0.315 

0.05 

11.0 

0.390 

0.05 

15.0 

0.475 

0.05 

22.5 

0.575 

0.06 

30.0 

0.665 

0.06 

37.5 

0.745 

0.07 

namely,  the3'  agree  fairly  well  with  the  line  measurements,  although  they  are 
slightlj’  higher.  R3T  measurements  from  a  different  wafer,  again,  agree  well 
with  those  of  R2T  and  the  line  data.  Still  the  ones  at  the  higher  frequencies 
are  more  near  the  approximate  formula  calculations  than  the  other  measured 
and  MDS  data. 

One  solid  conclusion  that  can  be  drawn  here,  at  least  for  relativeh' 
thick  gold  conductors  on  GaAs  substrates  which  are  typical  of  MMIC’s,  is  that 
some  of  the  approximate  formulas  available  will  overestimate  the  loss  by  almost 
a  factor  of  two  at  the  lower  frequencies.  Also,  MDS  will  slightly  underestimate 
the  loss  for  smooth  conductors  and  be  very  good  for  rough  conductors.  Of 
course,  GaAs  MMIC’s  will  typically  have  very  smooth  surfaces,  so  the  smooth 
behavior  of  MDS  is  of  more  interest  in  this  case. 


Figure  4.21  Comparison  of  the  measured  and  calculated 
parallel  losses 


Frequency  (GHz) 


Figure  4.22  Comparison  of  the  measured  and  calculated 
transverse  losses 


CHAPTER  V 


SOLUTION  OF  THE  INTEGRAL  EQUATION 
FOR  THE  CURRENT  DISTRIBUTION  IN  A  STRIP  CONDUCTOR 

5.1  Basic  integral  equation  solution  technique  [44] 

Given  a  homogeneous  eigenvalue  equation  of  the  form 

Lif)  =  XM{f)  (5.1) 

where  L  and  M  are  operators,  allowable  A's  are  eigenvalues,  and  the  corre¬ 
sponding  solutions  /  are  eigenfunctions,  the  method  of  moments  can  be  used 
to  reduce  this  functional  equation  to  a  matrix  equation.  An  inner  product  def¬ 
inition  if,g)  is  also  required;  this  inner  product  is  defined  to  be  commutative, 
associative,  and  positive  definite.  The  equation  (5.1)  often  arises  in  phj'’si- 
cal  problems;  the  eigenvalues  correspond  to  characteristic  parameters  (such  as 
propagation  constants)  of  these  problems. 

The  general  approach  to  solving  (5.1)  with  a  moment  method  is  to 
choose  a  set  of  basis  or  expansion  functions  /i ,  /2 ,  .fa ,  •  •  • ,  in  the  domain  of  L 
and  M,  and  expand  f  in  a  linear  combination  of  these  basis  functions: 

(5.2) 

n 

where  the  are  constants.  This  expansion  is  substituted  into  (5.1)  and  yields 

5]Q„£(/„)  =  A^a„M(.f„.)  (5.3) 

n  n 

We  take  the  inner  product  definition  mentioned  previously  and  choose  a  set  of 
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testing  or  weighting  functions  tf’i,  t{>2,  u’3, . . . .  in  the  range  of  L  and  M.  If  the 
’s  are  the  same  functions  as  the  /„, 's,  then  this  is  referred  to  as  Galerkin’s 
method.  The  inner  product  of  (5.3)  is  taken  with  each  Wm,  giving 

(Xfilu'fjj,  L  fn)  —  A  ^  ^  Ad fn)i  V??7..  (6.4) 

n  n 

This  can  be  rewritten  in  matrix  form  as 

[^mn][^ln]  ~  (5-^) 

The  equation  (5.5)  only  has  solutions  if 

det|/njj7  \  ~  O*  (5.6) 

The  determinant  is  a  polynomial  in  A  whose  roots  are  the  eigenvalues  of  the  ma¬ 
trix  equation  (5.5)  and  the  approximate  eigenvalues  of  the  functional  equation 

(5.1). 

The  operators,  of  course,  can  be  integral  or  integro-diffei'ential;  fur¬ 
thermore,  the  resulting  matrix  equation  may  be  a  nonlinear  eigenvalue  equation 
that  only  has  solutions  (which  are  often  difficult  to  find,  [47])  when 

det|F(A)|  =  0.  (5.7) 

This  is  the  case  for  the  integral  equation  presented  here. 

5.2  Problem  statement 

Figure  5.1  shows  the  arrangement  and  coordinate  system  to  which  the 
moment  method  will  be  applied.  An  integral  equation  which  is  applicable  to 
this  geometry  is  available  [45].  The  unknown  in  this  equation  is  the  current 
distribution  in  the  imperfectly  conducting  non-zero  thickness  strip  of  arbitrary 
cross  section  on  a  perfect  dielectric  above  a  perfect  ground  plane.  The  objective 
is  to  solve  for  the  attenuation  constant  a  for  this  structure.  To  simplify  matters, 
the  solution  will  first  be  carried  out  for  the  case  of  =  1  (air  dielectric).  Payt 
[46]  has  carried  out  similar  work  for  dielectric  waveguide. 


perfect  conductor 


Figure  5.1  Geometry  of  the  microsti'ip  problem 


5.3  Air  dielectric  case  (substrate  Cr  =  1) 


5.3.1  Discretization  and  application  of  Galerkin’s  method 

The  fundamental  mode  for  this  structure  with  air  dielectric  is  given 


E  =  E{x.y)e  '' 


H  =  Six. y)t- 


where  7  =  Q+jS  is  the  propagation  constant.  This  is  a  TEAI  mode  propagating 
in  the  direction  with  phase  constant  and  attenuation  constant  a.  An 
integral  equation  appropriate  for  this  structure  is  [45] 

r  =  (VV  •  +kl)  f  [Gj(x'.y')  +  GiJi(x',y')]  dx'  dy'  (5.9) 

Js 

where  E  =  ■  ,  and  (Tc  are  the  permittivitv  and  conductivity  of  the 

conductor,  e  =  €c—j^,  and  subscript  i  indicates  imaged  quantities.  The  kernel 
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is  a  modified  Bessel  function  of  the  second  kind,  such  that 

E  =  {VV-+kl)^  [  {ir-l)E{li')K,(u\p-p'\)dS' 

+  ( VV  ■  +1.-1)  ^  f  (e,  -  1)  £'■  (p',)  A%  (u\p  -  ;o;|)  dS’  (5.10) 

JSi 

where  ^  ~  p'  -  (a-',y'),  p'i  =  (-x',y'),  p'  €  5.and^'-  €  Si.  We  have 
set  u  equal  to  for  convenience. 

To  solve  the  integral  equation  (5.8)  we  grid  the  strip  conductor  of  Fig¬ 
ure  5.1  into  rectangular  cells  and  ex^Dand  E  and  Ei  in  a  finite  set  of  subdomain 
basis  functions  where  the  are  rectangular  pulses  which  are  nonzero  only  in 
their  respective  cells: 

M  M 

E^'^amfmiP);  Ei^  '^amjm(p)  (5.11) 

m  =  l  m=l 

and  let  the  testing  functions  also  be  fmiP)-  This  is  therefore  a  Galerkin  tech¬ 
nique.  Defining  the  symmetric  products 

(.f, 9)s  =  J^f-gdS]  if,4>)s  =  J  f<pdS  (5.12) 

and  adopting  the  notation 


(er-l)EKoi  )dS'={EAer-l)Ko{  ))s,,  (5.13) 

we  can  take  the  dot  product  of  f„{p)  with  the  equation  (5.3)  and  integrate  over 
S.  This  yields 


M  M 

53  ^m{LJn)s  =  ^  E  [(^r 

m=l  m=l 

M 

+  ^  E  {  (/-  •  +^^°] 


V  n  =  l,...,M. 


l).fm:.h.o(  )]5;)5| 

l)/m,A'o(  )]5,)5} 


(5.14) 
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For  the  rectangular  basis  functions  /„,  and  /„ ,  it  is  obvious  that 


{fm,fn)s  =  ^mn  ‘  area  of  rectangular  cell  =  Smrikp. 


(5.15) 


The  double  surface  integrals  in  (5.14)  are  now  converted  into  double 
contour  integrals.  For  brevity  we  will  examine  the  nonimaged  (subscript  m) 
portion  only.  We  begin  with 

■  |(VV-+A’2)^(eV-l).LA'o(t/lp-pT)  d5'|  dS  (5.16) 
and  let  K  =  (e^  —  l)/vo  {u\p  —  p'\)-  Then  we  can  write 
^  .fr,-|(VV-+A-2)^A'/„d5'| 

=  /  W.  /  luOfmdS'  ■f„ds+  [  kl  f  A-o()/„,d5'l  ■fr.dS 
Js  Js  J  Js  I  Js 

=  A  +  B.  (5.17) 


-4  =  ^V  J^V-KUdS'  ■f„dS  =  J^v\^£KiU-,l')dI'  ^f„dS  (5.18) 

by  direct  application  of  Stokes’  Theorem  to  the  integral  inside  the  brackets. 
Consider  next  the  vector  identitv 


V  ■{(f>D)  =  D-V(l)  +  (PV-D 


(5.19) 


and  let 


<l>^  £Kfr„-n'dr 


(5.20) 


Here  V  •  .D  =  V  •  /„  ,  but  /„  is  a  rectangular  pulse.  In  the  sense  of  generalized 
functions,  V  •  /„  is  a  delta,  function  at  the  left  side  of  the  pulse  and  a  negative 
delta  function  at  the  right  side;  as  this  is  inside  a  surface  integral,  we  can  (inside 
the  integral)  let  .D  •  V<^  =  V  •  i<f>I))  (since  the  integral  of  the  delta  functions  is 
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zero)  or 


.4  =  -n')  rf/'j  •  fn  dS 

=  ^V.  fn£K(.L-h’)dr  dS 


=  l^/n  K{fm  •  «')  dV  ■  n  dl 

=  ££R\L-r7Xfn-n)dl'  dl. 


(5.21) 


Now  for  the  B  portion  of  the  integral.  Since  K  is  a  Green's  function, 
it  satisfies  the  Helmholtz  equation  and  we  have 


( +  kl)K\  {ii.\p  -p'\)  =  (V2  +  kl  +  7')A-o  (u\p  -  p'\) 


^{V^  +  iXK.(u\p-p'\) 

=  -2t:S{p-p')\ 


(5.22) 


therefore 


u^Ko  {u\p-  p'\)  —  -  27: Sip  -  p')  -  V^Ko  (u\p  -  p'\)  =4> 

klKo  (u\p-p'\)  =  -  4  [+27r5(p  -  p')  +  V^A'o  iu\p-p’\)]  (5.2.3) 


Then 


B  =  4/f/ 

Js  Us 


i2nSip-p')  +  VU<)frr,dS'  ■fr.dS. 


(5.24) 


The  delta  function  portion  of  this  integral  is  straightforward: 


s  L-'s 


2t:  Sip  -  p')fmdS'  ■  fn  dS  =  27r/^  •  fn\p=p'  =  27r  S^n  (5.25) 


where  S^n  is  the  Kronecker  delta.  The  second  part  of  this  integral  is 


Is  [/4‘ 


K}fmdS'  -fndS 


IsMIs 


i^tK)-fmdS'  fndS 


=  (f  [  iVtK)  •  fm  dS'  ifn  •  h)dl.  (5.26) 
Jc  Us 
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We  make  use  of  the  vector  identity 

{VtK)  •  /to  =  ■  Kfm  -  KV  ■  fm  (5.27) 

as  before,  so  that  within  the  integral  V  •  fm  vanishes.  W'e  rewrite  the  integral 
as 


Vt  ■  Kfm  dS'  (/„  -  71.)  dl  - 


c  US 


Kfm  •  n'  dV  ifn  •  n)  dl 


c  Uc 


^  fcfc 
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Figure  5.2  Discretization  of  the  strip  into  a  mesh  of  rect¬ 
angular  cells 

5.3.2  Computer  implementation 

Figvure  5.2  shows  how  the  strip'conductor  was  gridded  up  into  a  mesh 
of  rectangular  cells.  The  cells  at  the  extreme  edges  of  the  conductor  are  square; 
this  allows  more  precise  representation  of  the  current  singularity  at  the  edges. 
Table  5.1  lists  the  basis  function  definitions,  while  Table  5.2  contains  the  dot 
product  evaluations  for  the  contours  shown  in  Figure  5.3.  These  dot  products 
are  the  (fm  •  n')  and  (/„  •  n)  terms  from  Equation  (5.31).  Restricting  the  cells 
to  being  rectangiilar  keeps  the  dot  product  terms  simple;  one  only  needs  to 
know  the  relative  directions  to  know  the  sign,  whereas  only  parallel  edges  give 
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Table  5.1  Basis  function  definitions 

•  There  are  P  rectangular,  nonuniform  cells 

•  P  =  ■  Ni 

•  There  are  M  =  ZP  basis  functions  (3  per  cell:  x,y,z  parts). 

•  The  rectangular  basis  functions  are  defined  as  follows: 


fzp-2  —  ^x- 

p  €  cell 

P 

fzp—1  ~  ^yi 

p  G  cell 

P 

fzp  =  d;, 

p  €  cell 

P 

fm  —  otherwise 


nonzero  contributions 

The  image  plane  ( Y  axis)  gives  the  following  symmetry  relations: 

^Xi\{-x',y’)  =  P^x\{x',y') 

^yi\(-x',y’)  ~P'y\{x’,y') 

^zi\(-x',y')  =  ~Ez\{x',y') 

These  relations  let  us  convert  the  Um-  terms  in  Equation  (5.31)  to  their  cor¬ 
responding  ttm  terms,  thus  rendering  the  unknown  current  (or  electric  field) 
distribution  vector  in  entirely. 

It  is  important  to  note  that  when  the  primed  and  unprimed  coordi¬ 
nates  are  a  common  cell,  edge,  or  vertex,  special  measures  must  be  taken.  In 
these  cases  the  arguments  of  the  modified  Bessel  functions  can  be  zero  over  part 
of  the  range  of  integration,  and  Ko  has  a  singularity  for  zero  argument.  In  fact, 
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fm  •  1^'  =  0  when  fm  or  n'  =  f 


•  ??,  =  0  when  /„  or  n  =  z 


When /to  =  X  =>  MOD(???  .  3)  =  1 

[  1,  \i  dV  = -dy' ^  dV^_i 
fm  •  =  <  —  1,  if  dV  =  dy'  =4»  dl'^_2 

[  0,  other  two  sides. 

I"  1,  if  dl  =  —dy  dh-4 
In  •  »  =  s  —1,  if  dl  =  dy  =>  dh-2 
I  0,  other  two  sides. 


When /to  =  y  ^  M0D(7?7  ,  3)  =  2 

'  1,  if  d/' =  da-' =4> 
fm  •  n'  —  <  —1,  if  dV  =  —dx'  ~dV^_i 
0,  other  two  sides. 

(  1,  if  dl  =  dx  d?2-3 
In  ■  n  =  ■!  —1?  if  dl  —  —dx  dli-i 
K  0,  other  two  sides. 


When/„ 


M0D(??7,3)  =  3 


fm  -n'  =  0 


f„  -11  =  0 


Note:  dl'  runs  over  the  entire  conductor,  whereas  dl  runs  over  the 


current  cell. 
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Figure  5.3  Nornenclature  for  evaluation  of  the  contour  in¬ 
tegrals  for  primed  and  unprimed  cells 

the  small  argument  form  of  this  function  is 

A’o(r)  -ln(^)  -  7  as  r  0  (5.32 

where  this  7  is  Euler’s  constant.  Any  numerical  integration  scheme  which  tries 
to  evaluate  this  function  with  a  zero  argument  will  fail;  however,  this  is  an 
integrable  singularity  [48,  49].  All  of  the  special  cases  (shared  edges  and  shared 
vertices  for  various  integration  directions)  were  performed,  and  the  results  are 
all  various  combinations  of  logarithmic  functions  of  various  cell  dimensions.  For 
the  case  in  which  an  edge  is  shared  and  integration  for  both  segments  is  in  the 


same  direction,  the  following  integTal  results: 
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In  |a’ —  a:-'!  d;r  do-' =  (6  —  a)^  ln|5  — a|  — ^  .  (5.33) 

For  the  case  of  a  shared  vertex  (and  colinear  segments)  and  integration  in  the 
same  direction,  this  integral  is  required: 

a’  —  a:' I  dxclx'  =  —  ^(6  —  o.)^ln|5  —  a\  —  •^(c  —  6)^ln|c  —  b\ 

+  ^(c  —  a)"ln|c  —  o|  +  |(c  —  b){b  —  a).  (5.34) 

When  one  of  the  segments  is  integrated  in  the  oi323osite  direction,  the  resulting 
integrals  are  simph*  the  negative  of  those  in  (5.33)  and  (5.34). 

The  code  for  the  air  dielectric  case  was  written  and  apparently  works, 
but  incompatible  compiler /operating  system  versions  at  my  computer  site  pre¬ 
vented  me  from  getting  results  in  time  to  include  in  this  paper.  Work  on  that 
item  is  proceeding  separately. 

5.4  Perfect  dielectric  case  (substrate  tr  >  1) 

There  was  not  enough  time  to  fully  address  this  case.  The  main  dif¬ 
ference  between  this  case  and  the  air  dielectric  case  is  that  the  Green's  function 
is  different.  It  now  involves  Sommerfeld  integrals  instead  of  contour  integrals; 
this  is  a  result  of  having  a  source  over  an  infinite  dielectric. 


CHAPTER  VI 


CONCLUSIONS 

This  dissertation  contains  the  results  of  work  carried  out  both  at  the 
University  of  Colorado  at  Boulder  and  at  the  Avionics  Laboratory  at  Wright- 
Patterson  Air  Force  Base  in  Dayton.  Ohio.  This  work  in'\'oh’’ed  fabricating  gold- 
on-GaAs  microstrip  lines  with  tightly  controlled  surface  roughness  and  cross 
sectional  area,  as  well  as  making  careful  on-wafer  measurements  of  these  lines 
attenuation  constants.  More  work  was  performed  on  the  solution  of  an  integral 
equation  which  can  be  used  to  calculate  loss;  however,  this  particular  equation 
is  not  usually  used  for  this  purpose  since  applying  the  method  of  moments 
yields  a  nonlinear  complex  eigen\-alue  problem.  The  complex  eigenvalues  are 
the  complex  propagation  constants  for  the  modes  which  may  exist  on  the  line. 
We  were  only  interested  in  the  dominant  mode. 

Much  work  was  expended  in  developing  and  characterizing  the  process 
used  to  fabricate  the  passive  MMIG's  for  this  project.  Many  of  the  techniques 
could  be  used  here  at  the  University  of  Colorado  with  minor  changes  in  the  way 
things  are  done  now.  Making  these  changes  and  documenting  the  procedures 
would,  I  believe,  greatly  improve  any  in-house  processing  of  GaAs  which  might 
be  carried  out  here.  I  hope  that  this  document  is  the  start  of  such  improvement. 

Much  work  was  also  expended  in  setting  up  the  measurement  system 
here  at  the  Universitj"  of  Colorado  and  in  getting  good,  useful  data  from  this 
system.  This  work  has  helped  provide  the  university  with  the  ability  to  make 
on- wafer  measurements  of  both  passive  and  biased  active  devices  up  to  49  GHz. 
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Again,  many  of  these  details  and  procedures  are  outlined  in  this  work,  and  it 
will  hopefully  be  a  useful  guide  to  those  who  come  later.  The  use  of  on- wafer 
TRL  calibration  and  the  standards  thereof  are  also  discussed  in  some  detail. 
The  data  collected  and  presented  show  that  this  system  can  make  good  mea¬ 
surements.  Some  work  on  computing  and  characterizing  the  effects  of  different 
errors  on  the  TRL  calibration  needs  to  be  done.  In  hindsight.  I  believe  that 
work  expended  on  the  integral  equation  solution  would  have  been  more  logicallj’ 
spent  on  a  thorough  error  analysis  of  the  effects  of  real-world  variations  on  TRL 
calibrations;  effects  such  as  via  hole  variation,  nonuniform  substrate  height,  and 
differences  in  probing  are  not  going  to  go  away  and  must  be  accounted  for. 

Finally,  work  was  started  on  the  integral  equation  solution.  Working 
code  was  written .  for  the  air  dielectric  case,  but  a  compiler  problem  on  my 
computer  system  prevented  running  the  code  enough  to  get  results  to  include 
in  this  dissertation. 
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APPENDIX  A 


INITIAL  WAFERS 

A.l  Non-via-liole  coplaiiar-to-microstrip  transition 

At  the  start  of  this  research  there  was  a  lack  of  processing  experience. 
Therefore  the  decision  was  made  to  use  a  coplanar  prolse  to  microstrip  transition 
(Figure  A.l )  that  does  not  require  via  holes  [19],  therefore  saving  one  mask  and 
two  processing  steps  (via  lithography  and  via  etch).  This  transition  is  a  radial 
stub  whose  dimensions  are  selected  so  that  its  impedance  to  (microstrip)  ground 
is  minimized  over  a  certain  frequency  range.  Of  course  the  ground  contacts  of 
the  probe  make  contact  with  this  stub,  and  the  center  conductor  of  the  coplanar 
waveguide  contacts  the  microstrip  line  proper. 

There  are  three  major  drawbacks,  though.  First,  the  structure  will 
not  have  a  constant  impedance.  Figures  A. 2  and  A. 3  show  the  theoretical 
impedance  of  the  transitions  used  in  the  first  two  wafer  designs.  Second,  these 
structures  consume  large  amounts  of  area  on  the  wafer.  Due  to  the  dependence 
of  dimensions  on  frequency,  lower  frequencies  require  prohibitively  large  stubs. 
Increasing  substrate  height  h  also  increases  size.  Third,  by  their  verj'  nature 
these  transitions  are  resonant  or  neax-resonant  structures.  If  one  is  trying  to 
measure  loss  and/or  measure  resonator  responses,  then  there  may  be  difficulties 
due  to  these  extra  resonances  in  the  measurement  path.  Calibration,  however, 
should  remove  concern  for  the  first  effect  and  remove  or  reduce  concern  for  the 
third  effect.  The  second  effect  is  a  physical  phenomenon  and  must  simply  be 
tolerated  as  an  engineering  tradeoff  for  not  having  to  do  via  hole  processing. 
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Figure  A.l  Coplauar  probe  to  microstrip  transition  •vi’ith- 
out  vias 


For  completeness  the  design  equations  as  derived  in  the  reference  [51- 
53]  axe  presented  below.  The  dimensions  refer  to  those  shown  in  Figure  A.l, 
and  X  is  the  reactance  of  the  stub. 
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Figure  A. 2  Impedance  of  non-via-hole  coplanar  transition 

for  first  wafer  design 
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Figure  A. 3  Impedance  of  non-\’’ia-hole  coplanar  transition 

for  second  wafer  design 
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In  the  calculations  for  these  wafers  ri  was  selected  on  basis  of  a  5  mil  probe  tip 
pitch. 


A. 2  First  wafer  designs 


Figures  A. 4  and  A. 5  are  views  of  the  first  two  wafer  designs.  Only  one 
of  each  wafer  was  successfully  fabricated.  Some  initial  measurements  were  made 


wafer  design 
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on  the  first,  and  the  second  was  broken  shortly  after  completion.  The  first  wafer 
is  a  100  fJ.  thick  design  just  like  the  newer  grooved  wafer.  Some  two-port  and 
one-port  gap-coupled  resonators  are  included,  and  their  resonant  frequencies 
were  selected  to  be  around  12  GHz  and  17  GHz.  There  were  long  lines  and 
calibration  standards.  This  wafer  was  used  in  setting  up  the  probe  station  at 
the  University  of  Colorado.  The  second  wafer  had  only  one  set  of  calibration 
standards  and  just  resonators.  The  noteworthy  thing  about  this  wafer  is  that  it 
was  designed  to  be  fabricated  on  unthinned  (20  mil)  GaAs,  therefore  eliminating 
another  fabrication  step.  Thus  these  transitions  are  bigger  and  spaced  farther 
apart.  This  wafer  was  simply  an  attempt  to  have  something  which  could  be 
easily  and  quickl}?^  fabricated,  which  seems  to  be  something  which  cannot  be 
asked  of  GaAs. 
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MEASUREMENT  SETUP 

This  appendix  presents  the  measurement  setup  at  the  University  of 
Colorado  at  Boulder,  and  discusses  how  to  use  it.  This  setup  is  very  versatile 
and  provides  the  following: 

1)  general  network  analyzer  functions, 

2)  optional  computer  control  of  network  analyzer, 

3)  on-wafer  probing  of  passive  devices,  and 

4)  on- wafer  probing  of  active  devices. 

Figure  B.  1  is  a  block  diagram  of  this  measurement  system. 

B.l  Power  Up  Sequence  and  Station  Assembly 

Ideally  the  851  OB  and  associated  equipment  would  never  be  turned 
off.  This  would  minimize  thermal  drift,  which  is  a  dominant  source  of  error. 
However,  the  equipment  usually  gets  turned  off  for  various  reasons,  so  the  first 
thing  to  do  when  making  measurements  is  to  turn  everything  on  and  allow  it  to 
warm  up  for  at  least  an  hour  -  preferably  several  hours.  The  power  up  sequence 
is  not  critical  except  for  one  thing:  the  851  OB  display  and  processor  unit  must 
be  turned  on  last;  otherwise  there  may  be  bus  contention.  One  useful  sequence 
is  to  turn  on  equipment  in  this  order: 

1)  Outlet  strip  A 

2)  Microscope  light  (setting  2  is  a  good  compromise  between  brightness 

and  heat  generation) 

3)  Chuck  control  power  supply 
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4)  S341B  source  (from  standby) 

5)  8514A  test  set 

6)  Outlet  strip  B 

7)  HP  9000/236  computer 

8)  851  OB  display  and  processor 

Finallj’’,  the  dual  3-1/2'’  disk  drive  can  be  connected  to  two  places.  If 
the  computer  will  be  controlling  the  measurements  (i.e.,  if  NIST's  DEEMBED 
program  is  being  used)  and  the  floppy  is  to  be  accessed  from  the  computer, 
then  its  cable  must  go  to  the  HPIB  connector  on  the  back  of  the  hard  drive.  If 
the  851  OB  is  being  used  and  flies  will  be  stored  and/or  recalled  directly  from 
the  front  panel,  then  the  cable  should  be  attached  to  the  connector  on  top  of 
the  plotter. 

Due  to  the  large  number  of  students  with  access  to  this  area  and  the 
tips’  fragility,  probe  tips  are  not  left  on  the  probe  station.  This  is  unfortunate, 
as  it  increases  chances  of  probe  tip  breakage  through  increased  handling.  It 
also  greatly  increases  wear  on  the  probe  tip  and  cable  connectors,  which  must 
be  put  on  and  taken  off  each  time. 

To  mount  the  probe  tips,  make  sure  the  probe  arms  are  backed  about 
halfway  back  from  the  stage  (using  the  x-axis  knobs  on  the  left  and  right  ends 
of  the  probe  arms).  This  gives  enough  room  to  turn  the  Allen  wrench  with¬ 
out  hitting  the  microscope  housing.  Now  mount  each  probe  tip  in  turn  using 
the  screws  and  an  Allen  wrench.  Be  careful  not  to  let  the  probe  tips  contact 
anything,  and  tighten  them  down  securely.  Now  attach  the  cables  to  the  probe 
tip  connectors.  The  cables  are  flexible  to  a  certain  degree.  Grasp  the  cable 
and  hold  it  so  that  the  threads  line  up  with  those  on  the  probe  tip  connector, 
and  tighten  it  down  finger-tight.  Using  the  3.5mm  torque  wrench  and  keeping 
the  cable  from  rotating  by  holding  a  finger  in  front  of  it,  tighten  the  connector. 
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Now  the  station  is  ready  for  sample  loading  and  probing.  When  the  measure¬ 
ment  session  is  concluded,  shut  down  the  systems  and  remove  the  probe  tips 
b}"  reversing  the  order  of  the  above  procedure. 

B.2  Mechanics  of  Probing 

Before  one  attempts  probing,  one  should  have  read  the  documenta¬ 
tion  contained  in  the  probe  station  user’s  manual  and  one  should  have  recei'^’ed 
hands-on  training  from  someone  with  prior  probing  experience.  The  documen¬ 
tation  contained  in  the  other  iDarts  of  this  manual  describe  the  operation  of  the 
probe  arm  micrometers  and  the  electronic  chuck  control  system. 

To  dispel  anj'  mj’ths  or  rumors  to  the  contrary,  the  z-axis  stage  or 
chuck  will  not  ‘hun  away’’  and  crush  the  sample  against  the  probe  tips.  If  there 
is  a  glitch  or  power  shutdown  and  restoration,  the  stage  will  go  through  its  auto¬ 
zero  routine.  This  routine  involves  the  chuck  going  to  its  zero  position,  which 
is  all  the  way  down.  As  it  does  its  calibration,  it  cycles  up  and  down  below  this 
zero  position.  Any  type  of  glitch  or  reset  will  cause  the  chuck  to  go  down-not 
up.  The  only  way  to  damage  the  sample  or  the  probes  is  operator  error — by 
manuallj’  running  the  stage  up  against  the  probe  tips.  The  HP  computer  does 
not  control  the  z-axis  stage! 

Probing  a  sample  requires  a  plan.  One  should  know  beforehand  what 
things  will  be  measured,  what  standards  will  be  used  for  calibration,  etc.  and 
have  a  sequence  in  mind.  This  minimizes  the  time  for  the  test  set  to  drift. 
The  normal  time  scale  after  which  drift  noticeably  affects  measurements  (es¬ 
pecially  using  on-wafer  TRL  and  higher  frequencies,  above  30  GHz)  is  about 
two  hours.  Since  placing  the  probes,  saving  data,  and  computing  coefficients 
(using  NIST’s  DEEMBED)  are  very  time-consuming,  it  is  mandatory  to  have 
a  probing  sequence  in  mind. 

When  ready  to  probe,  pull  the  stage  out  by  sliding  it  along  the  grease 
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plate.  Applj'  a  downward  as  well  as  a  horizontal  component  of  force;  if  an 
upward  component  is  applied  in  haste  or  carelessness,  the  grease  plates  could 
part  and  the  stage  could  fly  up  and  impact  the  probe  tips  (this  is,  of  course, 
operator  error!).  It  is  also  adviseable  to  move  the  stage  horizontally  slowh^ 
This,  like  all  facets  of  on-wafer  probing,  requires  practice.  After  the  stage  is 
out,  rest  the  sample  on  it  in  the  center  and  apply  vacuum.  Move  the  stage  back 
in. 

Using  the  microscope  on  about  SOX.  place  the  scope  in  the  center  and 
position  the  probe  tips  within  the  field  of  view.  This  is  somewhat  more  critical 
in  the  y-axis  (in-and-out)  so  that,  as  the  probes  are  moved  in  x,  they  remain  in 
view.  Now  focus  on  the  bottom  of  the  probe  tips.  The  reason  for  this  is  that 
when  the  sample  comes  into  focus  as  the  stage  is  raised  in  z,  it  will  be  in  contact 
when  it  is  in  focus.  Only  rarely  adjust  the  z-axis  micrometers!  Especially  leave 
the  LEFT  one  at  “0’‘.  If  the  probe  tip  bottoms  are  not  at  equal  heights,  adjust 
the  RIGHT  micrometer.  The  best  way  to  do  this  is  to  contact  a  noncritical 
area  of  the  sample. 

The  probe  tips  are  mounted  at  about  an  11-degree  angle  with  the 
horizontal.  After  the  tips  make  contact  and  the  stage  is  advanced  past  that 
stage  in  z,  the  tips  will  skid  on  the  surface  of  the  sample.  The  extra  height  jDast 
contact  is  called  overtravel,  and  due  to  the  1:5  angle  of  the  probes,  the  tips  will 
skid  5  units  on  the  sample  for  every  1  unit  of  overtravel.  Make  sure  to  use  the 
recommended  amount  of  overtravel  for  the  probe  tips  being  used. 

There  is  one  final  adjustment  to  be  made  before  measurements  can  be 
made,  and  it  must  be  checked  and  adjusted,  if  necessary,  before  every  contact  on 
the  sample.  This  adjustment  is  planarity.  The  only  way  to  effectively  tell  if  the 
probe  tips  are  planar  is  to  contact  the  sample  (in  an  unobtrusive  place)  NEAR 
where  the  contact  will  be  made.  If  only  two  marks  are  made  (or  one!)  then  the 
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probe  tip  must  be  rotated  with  the  supplied  Allen  wrenches  that  are  kept  on  the 
probe  arm  near  the  planarity  adjustment  screws.  Since  the  sample  will  rarely 
be  flat,  this  check  must  be  iDerformed  all  over  the  sample  if  accurate  results  are 
to  be  obtained  -  especially  with  TRL  on- wafer  calibration  of  microstrip  media. 

On-wafer  probing  is  a  skill  that  can  only  be  improved  with  practice. 
Everyone  develops  their  own  techniques  and  their  own  ways  of  “feeling”  what 
is  happening  with  the  equipment.  Care,  though,  is  always  required-care  and 
awareness  of  a  myriad  of  tiny  details.  Perhaps  this  appendix  will  eliminate  some 
of  the  tedium  and  time  of  flnding  and  learning  some  of  those  details. 
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CITIFILE  FORMAT  CONVERSION 

This  appendix  describes  how  to  port  data  from  the  HP  8510  network 
analj  zer  to  the  HP  workstations  and,  if  desired,  from  there  to  a  DOS  machine.  ^ 
The  conversion  of  the  resulting  files  to  a  format  useable  by  other  programs  or 
user-written  software  is  also  discussed. 

The  8510  currently  has  a  dual  disk  drive  attached  to  it;  disk  0  is  the 
one  to  which  the  8510  defaults  on  power-up.  This  drive  (and  the  ones  connected 
to  the  workstations)  only  read  and  write  730K  3-1/2"  floppy  disks. 

V  hen  you  want  to  save  data  tvhich  is  being  measured  and / or  displayed 
on  the  8510,  insert  a  disk  into  drive  0.  Press  the  TAPE  hardkey  at  the  bottom 
center  of  the  8510  display  unit.  There  are  two  softkeys  under  STORAGE  IS: 
they  are  TAPE  and  DISC.  If  the  DISC  selection  is  not  already  indicated,  press 
that  softket-.  (The  machine  will  power  up  with  storage  set  to  TAPE.)  Now 
press  the  STORE  softkey;  a  menu  of  choices  as  to  what  data  can  be  stored  will 
present  itself.  Usually  you  will  want  to  press  MORE  and  DATA  (under  the  DATA : 
heading.  The  measured  data  will  be  stored.  The  three  choices  under  the  DATA : 
heading  are  RAW,  DATA,  and  FORMATTED.  Raw  data  is  what  the  8510  measures 
before  calibration  is  applied;  data  under  the  DATA  heading  has  been  processed 
with  the  current  cal  set  (if  calibration  is  on).  Formatted  data  is  stored  in  the 
manner  in  which  it  is  displayed.  (Example:  if  j'^ou  were  viewing  Sn  data  in  log 

^  Step-by-step  commands  pertain  to  the  network  in  the  MIMICAD  Cen¬ 
ter  at  the  University  of  Colorado  at  Boulder. 
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magnitude  format,  it  will  be  stored  as  log  mag  Sn,  0.0)  Data  is  normally  saved 
in  real,  imaginary  format.  Formatted  data,  though,  is  just  what  is  displayed 
upon  the  S510's  screen.  The  8510  stores  these  files  in  a  format  called  “lif".  It 
also  capitalizes  the  file  names.  Since  UNIX  is  case-sensitive,  these  filenames 
must  be  tj'ped  in  upper  case  on  the  workstation.  Once  the  data  files  are  stored 
on  disk,  they  can  be  cojoied  to  your  workstation  account. 

There  are  two  HP  workstation  drives,  referred  to  as  /dev/rmdleft : 
and  / dev/ rmdright : .  These  drives  are  connected  to  kamikaze,  so  if  you  logged 
on  to  a  different  workstation,  you  must  rlogin  kamikaze.  Once  you  are  prop¬ 
erly  logged  in, 

lifcp  /dev/rmdleft : source-file  destination-file 
or 

lifcp  /dev/rmdright: source-file  destination-file 
will  copy  the  desired  file  to  your  account.  Other  basic  file  commands  are  made 
by  joining  “lif”  with  basic  UNIX  file  commands;  for  instance, 
lifls  /dev/rmdleft: 

will  list  the  files  on  the  disk  in  the  left  drive. 

UNIX  and  DOS  store  files  in  a  slightly  different  manner.  In  order  for 
other  programs  (such  as  editors  and  plot  programs)  to  work  correctly  with  these 
files,  a  conversion  must  be  performed  after  a  transfer  in  either  direction.  If  a 
file  is  to  be  copied  to  a  DOS  disk,  it  must  be  converted  with 
ux2dos  fileneone  >  newname. 

On  the  other  hand,  if  a  file  was  just  copied  from  a  DOS  disk  (or  uploaded  via 
modem),  this  command  must  be  executed: 
dos2ux  filename  >  newname. 

When  you  are  ready  to  copy  files  to  a  DOS  disk,  commands  similar  to 
lifcp  and  lifls  are  available.  The  command 
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doscp  source-file  /dev/rmdright : destination-file 
copies  a  file  from  the  workstation  to  a  DOS  disk  in  the  left  drive.  Remember, 
DOS  onh'  allows  eight  characters  in  a  filename  and  three  characters  in  the 
extension,  and  is  not  case-sensitive.  Likewise. 

dosls  /dev/rmdright: 

lists  the  files  on  a  DOS  disk  in  the  right  disk  drive. 

The  HP  8510  also  arranges  the  data  in  the  file  it  stores  in  a  manner 
referred  to  as  citifile  format.  This  refers  to  the  arrangement  of  information  in 
the  header  and  the  ordering  of  data  in  the  rest  of  the  file.  Figure  C.  1  is  an 
example  of  a  raw  data  file.  The  first  two  lines  are  always  the  same  for  this 
machine.  The  next  line  indicates  whether  data  is  raw,  measured,  or  formatted 
(NAME_RAW_DATA,  NAME-DATA,  or  NAME-FORMATTED,  respectively).  The  next  line 
gives  some  information  not  normally  needed.  In  the  next  four  lines  the  data  that 
is  actually  stored  is  indicated;  these  are  usually  S-parametefs  in  real,  imaginary 
(RI)  format.  The  middle  line  of  the  next  three  lines  always  contains  SEG,  a 
space,  start  freciuency  in  hertz,  a  space,  end  frequency  in  hertz,  a  space,  and 
the  number  of  points.  Each  S-parameter  is  then  stored  with  a  line  saying 
BEGIN,  then  the  S-parameters  themselves,  and  then  a  line  saying  END.  This  is  a 
file  format  which  is  of  little  use  to  most  other  programs. 

A  FORTRAN  program  called  convert. f,  compiled  to  an  executable 
convert,  is  available.  It  prompts  for  source  filename,  checks  to  see  if  the  source 
is  a  valid  data  file,  and  finds  out  how  many  S-parameters  are  stored.  The  user 
then  enters  the  destination  filename  and  indicates  whether  separate  files  or  one 
multi-column  file  are  desired.  Separate  files  will  have  appropriate  extensions 
(.Sll,  .321,  etc.).  The  one  combined  data  file  will  have  the  extension  .DAT, 
but  if  this  option  is  selected  and  only  one  S-parameter  is  stored,  it  will  still 
be  given  the  proper  .Smn  extension.  The  program  stores  converted  data  in 
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columns  with  frequency  in  GHz  in  the  first  column  and  converted  S-parameters 
in  the  following  columns,  in  one  of  a  choice  of  formats:  log  magnitude/phase 
(degrees),  linear  magnitude/phase  (degrees),  or  real/imaginary.  The  source 
code  is  available  and  easily  modified. 


Table  C.l  Example  HP  8510  raw  data  file  in  citifile  format 


CITIFILE  A. 01. 00 

#NA  VERSION  HP8510B.05.00 

NAME  RAWJDATA 

#NA  REGISTER  1 

VAR  FREQ  MAG  801 

DATA  SCl,l]  RI 

DATA  S[2,l]  RI 

DATA  S[l,2]  RI 

DATA  S[2,2]  RI 

SEG_LIST_BEGIN 

SEG  15000000000  18000000000  801 

SEG_LIST_END 

BEGIN 

4 . 97345E-1 , -4 . 86526E-1 
3 .40637E-1 ,-6 . 30218E-1 


5.16143E-1,-4.17022E-1 
3 . 29437E- 1,-5. 48584E- 1 
END 
BEGIN 

-2 . 54839E-4 , -0 . 03740E-4 
-4 . 33698E-4 , -4 . 10512E-4 


-3.15159E-4,-5.86032E-4 
-5 . 95986E-4 , -7 . 71 135E-4 
END 
BEGIN 

-2 . 36421E-4 , 0 .43705E-4 
0 . 17419E-4 , -2 . 93001E-4 


-1 . 68502E-4 , -1 . 14297E-3 

-0.85659E-4,-1.60181E-4 

END 

BEGIN 

-5.18096E-1,5.07598E-1 
-3 . 30230E-1 , 6 .30218E-1 


-3.36669E-1,5.16082E-1 

-1.83136E-1,6.15173E-1 

END 


